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ABSTBACT 


An investigation of the passive shock wave/boundary layer 
interaction control aiming at reducing the drag for conventional and 
supercritical airfoils at transonic Mach numbers is presented. A 3x 
lS.4-lnch Transonic Wind Tunnel was designed, constructed and calibra- 
ted to achieve this objective. Modifications were made in the Initial 
constant area test section to accomodate for the boundary layer growth 
along the tunnel walls. The boundary layer of the test section bottom 
wall was removed via a bleed system, so that the new boundary layer 
began at the airfoil leading-edge stagnation point. A variable porosity 
test section top wall was used to minimize the wall Interference. A 
manometer board and a Schlleren system were contructctd to measure the 
pressures and obtain Schlleren photographs of the flow field over the 
different airfoils In the test section. 

The passive drag control concept, consisting of a porous sur- 
face with a cavity beneath It, was Investigated with a 12-percent -thick 
circular arc and a 14-percent-thlck supercritical airfoil mounted on 
the test section bottom wall. The porous surface was positioned in the 
shock wave/boundary layer Interaction region, the flow circulating thr- 
ough the porous surface, from the downstream to the upstream of the 
terminating shock wave location, produced a lambda shock wave system 
and a pressure decrease In the downstream region minimizing the flow 
separation. The wake Impact pressure data showed an appreciably drag 
reduction with the porous surface at transonic speeds. To determine 
the optimum size of porosity and cavity, tunnel tests were conducted 


xix 


with different elrfoll porosities, csvltles and flow Mach numbers. A 
higher drag reduction was obtained by the 2.5 percent porosity and the 
l/A-lnch deep cavity. 
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PART 1 


INTRODUCTION 

The Increasing cost and future scarcity of fuel for transport 
and military aircraft create an urgent need to reduce fuel consumption 
and increase aircraft performance through refinements in aerodynamics* 
lighter materials* and more efficient Jet engines. Several methods of 
reducing fuel consumption and improving performance are either being 
planned or under study for existing and future aircraft. The reduction 
of wing drag is a logical method of achieving this goal* especially at 
the transonic cruise speeds of Jet aircraft. 

In this transonic flight regime, where complex nonlinear ef- 
fects tend to dominate the flow, a dramatic Increase in the drag takes 
place as the drag divergence Mach number is reached. The most common 
description of a transonic flow is when there is a supersonic “bubble" 
totally imbedded in a subsonic flow. The supersonic bubble may be ter- 
minated by a shock wave producing wave drag, or in certain circumstances* 
may return to subsonic conditions through an isentropic compression with 
no wave drag. This wave drag is associated with the increase in entropy 
across the shock wave. Soon after a shock wave appears in the flow* the 
drag will Increase rapidly with increasing free-stream Mach number lead- 
ing to the "drag rise Mach number," which is defined as the free-stream 
Mach number at which this rapid drag rise begins. One of the main objec- 
tives of designing a wing for transonic speeds is to obtain as high a 
"drag rise Mach number" as possible* subject to certain constraints. 

The obvious way to reduce the wave drag* at least in the 2- 
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dimensional case, la to use a supercritical shock>free airfoil section^ 
where there is no shock wave and consequently, no wave drag. However, 
these shock-free airfoils may have undesirable off-design characteris- 
tics, such as strong shock waves suddenly appearing tdien the Mach number 
is perturbed slightly from its design value. Thus, an important con- 
straint in the design of airfoils is that there should be good off- 
design behavior. 

In wing design with 3-dlmensional flow, the design can be 
altered by spanwise changes in the sweep and twist. This process, how- 
ever, complicates the design procedure. In principle, the supercritical 
airfoils are shaped to reduce the drag associated with energy losses 
caused by shock waves and flow separation, but still there is no generally 

applicable technique available which gives a net reduction of the total 

2 

drag at transonic speeds. 

For the thick, supercritical airfoils for the proposed flying 
load-carrying wing with airfoil thickness of 20 to 25 percent of the 
chord, the increase in the drag at transonic flight Mach numbers is more 
severe than for the thinner (approximately 10-percent-thlck) supercriti- 
cal airfoil wing. The drag of such a thick wing is reasonably low at 
the design lift coefficient and Mach number, but the drag increases at 
higher lift coefficients and Mach number because of the strong shock 
waves which terminate the imbedded supersonic region. 

In the case of fighter aircraft, the performance envelope for 
maneuvering is limited by the onset of- buffet at high subsonic Mach num- 
bers, which is caused by severe unsteady shock wave /boundary layer inter- 
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action and Its effacta on the shock wave location. The tetnlnatlon shock 
waves for the Imbedded supersonic regions are very strong at high angles 
of attack and cause the boundary layer to separata from the wing. 

To control the drag Increase due to the shock wave/boundary 
layer Interaction for conventional airfoils and for thin and thick super- 
critical airfoils at truisonlc Mach numbers, a basic research program on 
the passive shock wave/boundary layer control for drag reduction was 
suggested by Mr. Dennis Bushnell and Or. Richard Whitcomb at the NASA 
Langley Research Center. The concept of the passive drag reduction con- 
sists of having a porous surface with a cavity underneath at the shock 
wave location. The high pressure downstream of the shock wave will force 
some of the boundary layer flow Into the cavity and out ^ead of the 
shock wave. By this susthod, the boundary layer will chicken ahead of 
the shock wave and send compression waves Into the supersonic region, 
thereby decreasing the Mach number for the normal shock wave. By this 
shock wave/boundary layer interaction process, the Increase of entropy 
across the shock wave will be lower and the boundary layer flow separa- 
tion will be minimized. Both of these effects tend to decrease the drag 
at transonic speeds. 

In recent experiments on the supersonic Jet noise reduction, 

3 

Maestrello used e porous centerbody plug nozzle Jet with a porosity of 
about 2 percent to eliminate the shock waves for transonic and supersonic 
Jet exhaust velocities. This method provided a shock free flow without 
the "shock associated noise" and "screechli^" and significantly reduced 
the Jet noise. These results were encouraging and Indicated the possl- 
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blllcy of mlalnlzlng th« shock wavs lossas ovar the airfoils at transonic 
Mach numbers by placing a porous surface with a cavity underneath the 
shock wave location. 

The effects of surface humping and chordwlse slots on the shock 
waves over supercritical airfoils at transonic Mach numbers were Investi- 
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gated by Lee and Yoshihara. The purpose of the Investigation was to 
decrease the adverse effects of shock wave/boundary layer interaction on 
severe buffeting at transonic maneuvering flight conditions. Prelimi- 
nary results did Indicate a change In the strength and location of the 
shock waves, with a slotted configuration showing greater effects. Air- 
foil static pressure distributions and Schlleren photographs were ob- 
tained at transonic Mach numbers and different angles of attack, but no 
drag measurements were taken for these two concepts. 

Ocher Investigators have studied the effects of Injection and 
suction on the control of turbulent boundary layars through transverse 
and longitudinal slots, holes, and porous plates at subsonic Mach numbers 
over airfoils. Wilkinson^ investigated the effects of suction on the 
Incompressible turbulent boundary layer through closely spaced stream- 
wise slots. Uls drag and mean boundary layer velocity profile measure- 
ments showed chat the slotted surface had nominally the same suction 
characteristics as the porous surface. A comprehensive list of refer- 
ences on the effects of suction and Injection on the turbulent boundary 
layer is presented In Ref. 5. 

Another technique of reducing the drag at subsonic Mach num- 
bers is laminar boundary layer control, which keeps the flow laminar. 


Thl« can be achlavad by uaing dlifarcnt approaches. Thasa include suc- 
tion, favorable pressure gradient, wall cooling for gaaea, wall hastily 
for liquids, MHD forces, and compliant walls. Among these approaches, 
the moat promislug one is laminar flow control by the uae of suction, and 
suction through porous atrlpa seems to be moat appropriate.^ 

Recent advances in computational aerodynamics have made it pos- 
sible to compute many nonlinear problems, such as the special solution 
for shock wave/boundary layer interaction flow over airfoils by Nandanan 
at al.^ However, this analysis did not include the separation of the 
boundary layer and the experiments showed the need for the correct treat- 
ment of the shock wave/boundary layer interaction in transonic flow. The 
strong viscous interaction between the shock wave and the boundary layer 
in the transonic region of nixed subsonic and st^rsonic flow is still 
not adaqiiately modeled. Despite some pioneering efforts, the develop- 
ment of methods for calculating the separated flow is still in its in- 
fancy. For this, we can refer to one of the techniques in Ref. 8, which 
presents the calculation of the interaction of the inviscid external flow 
and the separated boundary layer. The method has been applied to the 
flow field over a symmetric, biconvex airfoil at zero angle of attack, 
but has yet to be verified for the lifting airfoils. 

Investigations of ilie shock wave/boundary layer interaction 
over supercritical airfoils to Improve the theoretical analysis capability 
for airfoils experiencing strong Invlscld-viscous flow interactions are 
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being conducted at the Lockheed-Georgia Conq>any and at the NASA Ames 
Research Center}^ Because of the trend toward designing t h ic k e r and more 
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highly loaded airfoils for transport aircraft cruising at transonic 
spaads. it has bacome Increasingly difficult to predict the parforaanca 
of these airfoils. 

Because of the lack of experimental data and other urgent needs 
for decreasing aircraft fuel consumption and Improving aerodynamic effi- 
clency» the present basic concept of drag reduction has bem investigated 
mostly through experimental studies. The purpose of the research was to 
Investigate the possibility of drag reduction for supercritical airfoils 
by applying the concept of the passive shock wave/boundary layer control. 
The experimental Investigation was established to satisfy specific objec- 
tives idilch have provided basic Information about advanced supercritical 
airfoil wing design and the passive drag reduction concept over transonic 
airfoils. These objectives are as follows: (1) to define the flow field 

for the shock wave/boundary layer Interaction over the bottom surface of 
the transonic wind tunnel* which simulates the supercritical airfoil pres- 
sure distribution with a solid surface; (2) to determine the surface 
pressure distribution and the total drag from the impact surveys with a 
solid surface; (3) to study the effects of the surface porosity at the 
shock wave/boundary layer interaction region on the flow field and the 
total drag from the impact pressure surveys; (4) to Investigate the 
effects of the size of the cavity located below the porous surface on 
the flow field and the drag; (5) to determine the total drag from wake 
impact surveys and the flow field at transonic Mach nia^rs; (6) to 
study the effect of humidity on the dr^ reduction and the pressure dis- 
tribution. Schlleren photograpus for the optical Information were ob- 
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talned to define the shock waves with and without the porous surface for 
the passive drag control concept. 

To conduct this research, a blowdown transonic wind tunnel 
with an atmospheric intake was constructed using the existing dryer, 
vacuum tank and vacuum pumps used for the Mach 3 Wind Tunnel at Rensselaer 
Polytechnic Institute. A pressure measuring system, optical Instrumen- 
tation, and different models were designed, constructed and used to cali- 
brate the flow In the test section. A preliminary mahogany test section 
with constant area was used and modifications were made to obtain nearly 
uniform flow over the region where the airfoil model was located. Both 
Schlleren photographs and pressure distributions over a 10-percent- thick 
double wedge and 12-percent-thlck biconvex: airfoils were obtained with 
the models mounted in the center of the test section. Similar data were 
obtained with a 3-lnch chord circular arc airfoil mounted on the bottom 
surface of the test section. The details of these preliminary experi- 
mental results are presented in Ref. 11. An aluminum test section was 
later designed and constructed to investigate the concept of the passive 
drag control. The design of this test section was based on the knowledge 
gained from the Initial mahogany test section, and modifications were 
made to overcome the observed shortcomings. A 4-lnch chord circular arc 
airfoil and a NASA 14-percent-thlck supercritical airfoil (Ref. 12), were 
investigated In this new aluminum test section with an adjustable porous 
top wall to study the effects of the porous surface and the cavity size 
on the shock wave/boundary layer interaction. Schlleren photographs, 
pressure data and Impact pressure wake surveys were obtained at subsonic 
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and transonic speeds. 


PART 2 


EXPERIMENTAL FACILITY 


2.1 Transonic Wind Tunnel 

A 3>lnch by 15.4 inch transonic wind tiumel was constructed 
parallel to the existing 4-lnch by 4-inch Mach 3 Wind Tunnel to investi- 
gate the possibility of reducing the drag of the supercritical airfoils 
by the application of the concept of the passive shock wave/boundary 
layer control. The various components of the Transonic Wind Tunnel^^ 
are discussed below and Figs. 1 and 2 show the overall arrangement of the 
tunnel. 


2.1.1 Entrance Section and Settling Chamber 


A rectangular cross section settling chamber with a height of 
33.9 inches, a width of 20 Inches, and a length of 48 Inches was con- 
structed from 0.75 inch thick birch plywood. To smooth the entering air 
flow, a bellmouth was attached to the entrance of the settling chamber 
and to the wall of the dryer as shown in Fig. 1. 

A honeycomb with a depth of 6 inches and hexagonal flow pas- 
sages of approximately 0.25 inch was placed at the entrance of the set- 
tling chamber as shown in Fig. 1. Fine screens were placed in front and 
back of the honeycomb section to break up the large eddies into small 
sizes and to make the flow uniform across the settling chamber. The 
length of the settling chamber was selected to permit the decay of the 
small scale turbulent eddies before the flow was accelerated in the con- 
traction section. After the Initial tests, another layer of birch ply- 








ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 


a Transonic Wind Tunnel 


Mach 3 and Transonic Wind Tunnels 


cQTiJil 




tpi^ 1 





1 B ® 




I^Jd 












wood was added to the settling chamber walls to strengthen the rigidity 
of the section. 

2.1.2 Contraction Section 

The dimensions of the entrance to the contraction section were 
20 inches wide and 33.9 inches high, and at the exit the corresponding 
dimensions were 3 Inches and 15.4 Inches as shown in Fig. 1. A gradual 
contraction shape was selected towards the exit of the section to achieve 
uniform flow across the cross section. Multiple layers of thin aircraft 
plywood, glued together with epoxy, were used to obtain the desired con- 
tours for the contraction section. Longitudinal ribs were used to stiffen 
the plywood pieces as shown in Figs. 1 and 2, and later on epoxy filler 
was added on the sides of the contraction section exit to make it more 
rigid. The inside walls were carefully sanded and varnished to produce 
smooth surfaces for the air to accelerate from low speed to transonic 
speed at the end of the contraction section. 

2.1.3 Test Section 

The dimensions of the original test section were 3 inches wide 
and 15.4 inches high with a length of 36.5 inches as shown in Fig. 3a. 

For the initial calibration tests, solid mahogany pieces were used for 
the top and bottom surfaces so as to simplify the changes for these sur- 
faces. The test section had been modified to accommodate for the bound- 
ary layer growth in the test section by opening the top and bottom walls 
at the exit by 0.7 inch and the side walls by 0.15 inch as shown in Fig. 
3b. The initial 20 inches of the bottom and top walls were parallel to 
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the tunnel axis. For the initial flow field investigation, airfoils 
were placed in the middle of the flat bottom surface. Later, an alu- 
minum porous wall test section was designed and constructed, cf . Fig. 4. 
Its design was based on the knowledge gained from the initial mahogany 
test section. 

Although blockage effects are produced in most two-dimensional 
wind tunnel testings, the full impact of their presence is very evident 
in the transonic regime. VIhen the free stream Mach number approaches 
unity, a slight change in area ratio produces a large change in the Mach 
number. This happens when a model is placed in a test section operating 
in the transonic regime. 

For several years, slotted or perforated test section walls 

have been proposed to overcome these difficulties.^^ In our aluminum 

test section, attempts at reducing or eliminating the model blockage 

effects at transonic speeds have been made through the design and con- 
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structlon of an adjustable porous top wall test section. 

In the middle of the flat bottom surface of the test section, 
a 4-lnch chord circular arc and a 14-percent-thick supercritical air- 
foil with the porous surface and cavity located below the shock wave/ 
boundary layer interaction region, were placed and investigated. 

Thick Incite side walls were used to permit the observation of 
the shock waves and the flow field over the top surface of the circular 
and supercritical airfoils. Aluminum plates on both side walls were 
used to minimize the wall deflection. A circular cut was made in the 
aluminum plates for the optical information. With this arrangement, it 



g.A Side vio\: of iiluininun test section 


was possible Co study Che shock wave/boundary layer loteracclon phenomena 
with and without the passive porous surface for the circular arc and 
supercritical airfoils over a range of free scream Mach numbers. 

The test section boundary layer on the bottom wall was reroved 
via a separate bleed system, so chat the new boundary layer would begin 
at Che airfoil leading-edge stagnation point. The boundary layer removal 
system consisted of a 3-inch transverse slot located at the model leading 
edge. A 2-lnch diameter pipe connected the bleed slot to the vacuum 
system downstream of the wind tunnel. The removal of Che boundary layer, 
Fig. 8b, was established by opening a valve separating the vacuum Inside 
the pipe from the air flow In the test section at the start of the test. 

Two adjustable wedges, one on the bottom surface of the test 
section and the ocher on the top, were used to vary the free stream Mach 
number from low subsonic to transonic speeds as shown In Fig. 3b. This 
technique of controlling the Mach number In the test section has been 
effective. 


2.1.4 Diffuser Section 

The diffuser section consisted of two SO-lnch long diverging 
ducts as shown In Figs. 1 and 2. The original cross section at the en- 
trance was 3 inches by 15.4 Inches and these dimensions were Increased 
to 3.3 inches and 16.8 Inches as a result of the modifications In the 
test section. Fig. 3b. The dimensions at the exit of the diffuser were 
15.25 Inches by 15.25 inches. Two layers of 0.75-lnch chick birch ply- 
wood were used to construct the diffuser sections In order to withstand 
the pressure differential across the diffuser surface and to minimize 
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the noise from the transition shock waves. Between the two diffuser 
sections a manual control valve was installed as shown in Figs. 1 and 
2 to control the free stream Mach nuutber. However, a better control 
of the Mach number was obtained with the two adjustable wedges at the 
exit of the test section as shown in Fig. 3b. The test results with 
transonic flow in the test section have Indicated satisfactory perform 
mance of the diffuser to meet the design objectives. 

|if A plywood transition piece was constructed to connect the 

square diffuser exit cross section to a circular 14.5-inch diameter cross 
^ section ahead of the quick acting 16-inch dlasieter pneumatic valve as 

shown in Figs. 1 and 2. The transition from the square to circular cross 
section was accomplished by an octagonal section. 

2.1.5 Quick Acting Pneumatic Valve 

To start the flow through the transonic wind tunnel, a 16-lnch 
diameter quick acting pneumatic valve was mounted between the end of the 
diffuser section and the flange for the 16-lnch diameter vacuum pipe as 
^ shown in Figs. 1 and 2. A large 6-inch diameter pneumatic cylinder was 

' attached to the valve for quick opening. House air supply at a pressure 

of 120 psl was connected to the 6-inch diameter cylinder through a 110- 
volt solenoid valve. This valve was actuated by a key sicunted next to 
I the manometer board, cf. Figs. 1 and 2, to open and close the 16-lnch 

diameter valve rapidly. Tests had confirmed the fast opening and closing 

r 

I 

of the valve. A timing circuit had been installed to control the valve 

opening and closing as well as to close the valves on the manometer board 

I 

for measuring the pressure after a desired flow duration. In the closed 







18 


position with 120 psla prossure on th« control cylinder, the large valve 
sealed with practically no leakage to the vacuum tank. 

2.1.6 Vacuum System 

Six vacuum pumps with 5>hp electric motors were used to evacu- 
ate the large vacuum tank located adjacent to the laboratory. All of 
these pumps could be used simultaneously to evacusLe the tank, and could 
be left operating during a teat. It was possible to evacuate the tank 
to 29.5 inches of mercury pressure difference in a short time. 

The Mach 3 Wind Tunnel with a square test section of 4 Inches, 
located parallel to the Transonic Tunnel, as shown in Fig. 2b, utilized 
the same vacuum system. The discharge from both of these tunnels was 
connected to the 16-inch diai&eter pipe. 

2.1.7 Air Dryer 

The Transonic Tunnel was attached to a large cross section 
silica gel dryer as shown in Figs. 1 and 2. A sliding door was opened 
to permit the outside air to go through the dryer before entering the 
Transonic or the Mach 3 Wind Tunnels. Heating elements were installed 
to heat the recirculating air through the silica gel bed to 300°F for 
reactivation. After the silica gel bed was reactivated, a water-cooled 
heat exchanger was used to bring the bed temperature to ambient conditions. 

2.2 Instrumentation 

2.2.1 Pressure Apparatus 


A mercury manometer board was constructed with 20 U- tubes to 
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MMur* ch« scatic prcMurM «t various locations along cha Transonic 
Wind Tunnal and over the airfoil nodaxa. The same nanonatars ware used 
CO meaaura the impact prasaura from the rake for the wake survey and cha 
total pressure from an impact prasaura probe installed in cha settling 
chamber ahead of the contraction section. The free set earn Mach number 
was obtained by taking the average of the vertical static pressures 
meaaurad on the side of the entrance of the test section ahead of the 
airfoil modal location. Ball valves ware used with cha 0- tubas for the 
Initial calibration of the wind tunnel; but the valves were later re-> 
plead by an electrically timed solenoid valve system, which is discussed 
later. The valves were opened before the flow was started and ware 
closed by the timer after the steady pressure was established. Plastic 
cubing was used to connect the static pressure orifices in the tunnel 
and over the isodels to the valves attached to the U-tubes. 

2.2*2 Boundary Laver Survey Apparatus 

Ine boundary layer thickness along the side wall and the bottom 
surface of the initial test section was measured with a small Impact 
pressure probe. The diameter of the stainless steel probe was 0.063 
inch with an inside diameter of 0.0365 inch. This probe was attached 
CO a micrometer in order to accurately determine the distance of the 
impact pressure probe from the tunnel surface. The boundary layer thick- 
ness was measured at three axial locations for the side wall and the 
bottom surface. 

2.2.3 Wake Survey Apparatus 


A single Impact pressure probe was initially used to survey 
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th« imkc downscrcam of th« oirfoll modal. Bocauao of Its bonding In the 
flow when the flow woo ootabllshed, it was replaced by a chick rake with 
ten impact preeaure probes. The average probe spacing was 0.181 Inch 
and the average length from the probe centerline was 0.5 Inch. The diem* 
eter of the stalnlcas steel probe was .039 inch. The rake centerline was 
located 1.75 Inches downstresm of cha model trailing edge. Later the 
thick r^ce was replaced by a thin one to minimize the blockage effect. 

The thin rake had 8 Impact pressure probes with 0.031 Inch Inside and 
0.040 Inch outside diameters. The average probe spacing was 0.127 Inch. 
The location on the bottom wall test section was the same as for the 
thick rake. The survey of the wake was done through a height of 1.70 
Inches from the bottom surface. Thi) rake caps were connected to the 
mercury manometer board for measuring Che Impact pressures within the 
wake. 


2.2.4 Optical System 

The optical system was a conventional single-pass Schlleren 
system consisting of a focussed zirconium light source, two 7.5 Inch 
diameter parabolic mirrors with a focal length Of 5.0 feet, two 9.5 Inch 
diameter flat mirrors, an adjustable knife edge at the focal point, and 
a film plate holder. A schematic of this optical system Is shown In 
Fig. 5. With this system, it was possible to obtain both Schlleren and 
shadowgraph photographs of the shock wave and the flow field over airfoil 
models in the test section. The shadowgraph photographs were taken with 
the focal plane shutter camera placed next to the test section window. 

In this report, mainly the Schlleren photographs will be presented. 
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2.2.5 Solenoid Valves and Timer System 

Twenty small electrical solenoid valves were used Instead of 
the previous ball valves at the manometer board as shown In Fig. 2. They 
were normally opened before the flow was started and closed automatically 
through the use of an adjustable electrical delay timer. At the sane 
time the valves were closed, the timer activated the camera shutter for 
the Schlleren system. 

2.2.6 Temperature Gage 

A temperature probe was Installed In the settling chamber to 
measure the total temperature of the flow. The probe was connected to 
an electronic digital display system which read the temperature and dis- 
played It on a small screen. The temperature gage was connected to the 
timer system. The temperature was taken at the same time when the valves 
for the pressure were closed and the camera shutter activated. By this 
system, the pressures, temperature and the Schlleren photographs were 
taken at the same test time. 

2.2.7 Humidity Gage 

A humidity gage was installed In the wall of the dryer chamber. 
An opening made on the dryer wall allowed the air to flow Inside the gage 
for the readings. The relative humidity was in the Ideal range for most 
of the tests. After the dryer was reactivated, the gage Indicated a low 
relative humidity. 

2.3 Design and Construction of Airfoil Models 


Different airfoil models were designed and constructed to meet 
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Che requirements for the calibration of the Transonic Wind Tunnel and 
the investigation of the drag reduction by the concept of the passive 
shock wave/boundary layer control. 

2.3.1 Dobule Wedge Airfoil 

A 3-lnch chord, symmetrical double wedge airfoil of 10-percent 

thickness. Fig. 6, which spanned the test section of the 3- Inch by 15.4- 

Inch Transonic Wind Tunnel, was designed and constructed. This model was 

chosen for the wind tunnel test section calibration because experimental 

14 

data for the same airfoil tested In similar conditions were available. 

The airfoil model was mounted In the middle of the test section. The 
pressures over the model and the Schlleren photographs of the shock waves 
were taken for a range of free stream Mach numbers. 

2.3.2 Biconvex Airfoil 

A 3-inch chord biconvex airfoil of 12-percent thickness. Fig. | 

6, was designed next and constructed. The model was tested in the middle I 

of the test section. The biconvex airfoil was chosen basically because 
It best represented the conventional wing airfoil and also because experl- \ 

mental data were available for comparison, in Refs. 15 and 16. The pres- 
sures over the model and the Schlleren photographs were obtained for 
transonic free stream Mach numbers. 

2.3.3 3-Inch Circular Arc Airfoil 

A 3- inch chord circular arc airfoil of 12-percent thickness. 

Fig. 6, was designed, constructed and mounted this time on the bottom 

surface of the test section to investigate the pressure and the shock | 

’ 

: 




10-percent-thick double wedge airfoil 
3-lnch chord 


12-percent-thick biconvex airfoil 
3-inch chord 



12-percent-thick convex airfoil 
3-lnch chord 


12-percent-thick circular arc airfoil 
4-inch chord 



14-percent-thlck supercritical airfoil 
4-inch chord 


Fig. 6 


Airfoils 
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wave and to slmulace Che shock wave/boundary layer InCeracclon over Che 
supercrlclcal airfoil. Similar daCa, as for Che previous models, were 
caken. 


2.3.4 4-Inch Circular Arc Airfoil 

A 4-inch chord circular arc airfoil of 12-percenc Chlckness, 
Figs. 6 and 7a, was consCrucCed and cesCed on Che boccom surface of Che 
Cesc secclon. The model pressure dlscrlbudon and Che Schlleren phoco- 
graphs of Che shock waves over che model were obcalned. Based on this 
Informaclon, a porous surface wich a cavlcy undemeach ic was poslcioned 
ac Che locaclon of Che shock wave/boundary layer inceracclon. Boch Che 
porous surface and che cavlcy were used co invescigace Che drag reduction 
by che concepC of che passive shock wave/boundary layer conCrol. The 
porous surface conslsced of 12 rows of 38 holes each. The .021-inch 
dlameCer holes spanned Che airfoil model and Che porosicy exceeded from 
60 CO 90 percenc of che chord. The change of porosicy was obcalned by 
varying che porous surface lengch chordwise and by making Che holes 
larger. Based on Che model surface area, che porosicy could vary from 
zero CO 2.17 percenc. Three differenC values of porosicy were invescl- 
gaced: 2%, 1.3% and 1% porosicy. The no-poroslCy case was obCained by 

plugging Che holes. 

Holes of 0.021 inch dlameCer were drilled perpendicular Co Che 
model surface, and pressure cubes were connecced Co Chese holes. Ten 
scaclc pressures were discrlbuced along Che model cenCerllne wich an 
average spacing of 0.37 inch. The airfoil model was carefully sanded 
Co provide an aerodynamically smooch surface. The Inlcial excenc of 
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Fig. 7a 12-percent-thick circular arc airfoil 
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Che cavicy under the porous surface was from 40 to 94 percent of che 
chord and spanned the airfoil model. Ics dimensions were 2.125 inches 
long, 3 Inches wide and 2 inches deep. The length of the cavicy was 
later reduced to the porous surface size with a variable depth going 
from 2 inches to 0.25 inch, as shown in Fig. 8a. Investigations were 
conducted to study the cavity depth effect on che drag reduction and 
the flow field over che model. 

2.3.5 Supercritical Airfoil 

The basic concepts of supercritical airfoil technology were 

developed before adequate theoretical design or analysis codes were 

available. They were developed through an experimentally iterative process 

by evaluating che experimental pressure distributions at design and off- 

design conditions and physically altering the airfoil profile to yield 

che best drag characteristics over a range of Mach numbers. This is 

what happened at NASA Langley Research Center over the past several years, 

where "phas I airfoils" were developed and criteria were recognized which 

12 

provided guidelines for the design of supercritical airfoils. Based 

on these criteria, a 10-percenC-thlck^^ and a 14-percent- thick supercrltl- 
12 18 

cal airfoil * were designed. The design lift coefficient was 0.7 for 
both airfoils. An iterative design process was used which consisted of 
altering the airfoil coordinates until the viscous airfoil analysis pro- 
gram, developed in Ref. 19, indicated that the design criteria for a 
shock-free flow over the upper surface had been satisfied. 

For aircraft flying Just below the speed of sound, a wing is 
at its "critical" Mach number when the air flowing over it accelerates 
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to the speed of sound because of the Bernoulli effect. The air then ] 

slows down, producing a shock wave chat separates the boundary layer from 

, 

Che wing. Increasing drag and decreaslx^ lift. Supercritical wings allow 

an aircraft to fly faster before this occurs. The large radius at the 

front of the wing speeds up the air to supersonic velocity sooner, and - 

the flatter top surface sustains the high speed longer over the wing with 

a more gradual speed reduction. The drag<>producing shock wave Is moved i 

further back along the wing and Is Inherently weaker, as shown In Fig. 10. 

18 

A lA-percent'-thlck NASA supercritical airfoil. Fig. 11, was 
used In the present work to Investigate the drag reduction by the concept 
of Che shock wave/boundary layer control. The construction of the model 
was based on 200 points derived from the coordinates for the upper surface 
presented In Table 1. A computer program was used to define the model 
profile. The program was then used on a computerized milling machine 
to obtain a 14-percent thick aluminum supercritical airfoil. The model 
had a chord of 4 inches and a span of 3 inches. Fig. 7b. The airfoil 
was carefully sanded to provide a smooth model surface. Orifices of < 

0.021-lnch diameter were drilled perpendicular to the airfoil surface to 
measure the model pressure distribution. Seventeen static pressure Caps 

* 

were spaced along the centerline chord. The average spacing was 0.25 ; 

Inch; however, in the leading edge region, the spacing decreased to 0.18 

Inch. 

For the preliminary tests, the model was Investigated with a 


smooth solid surface. Model pressure distribution, impact pressure wake 
survey and Schlleren photographs were obtained. Later, a porous surface 
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was Introduced by drilling holes of 0.021- Inch diameter on the model 
surface. It extended from 56 to 83 percent of the chord. The choice of 
the porous surface location was based on the Schlleren photographs and 
the pressure data obtained from the Investigation of the model with a 
solid surface. The porous surface location was chosen such that the 
shock wave/boundary layer Interaction region was In the middle of the 
porous surface. The porosity » based on the model surface, varied from 
0 to 2.5 percent. The variation In the porosity was obtained by drilling 
larger holes or plugging some of them. The solid model surface was ob- 
tained by sealing all the holes. The porous surface spanned the airfoil 
model. Two different porosities of 0.94% and 2.5% were selected to In- 
vestigate the drag reduction and the optimum porosity size. The cavity 
used for the boundary layer flow circulation through the porous surface 
had the same length and width as the porous surface and Its depth could 
vary from 0 to 0.75 Inch. Two cavity depths 0.75 and 0.25 Inch were 
selected to study the cavity depth effect on the shock wave/boundary 
layer Interaction and to investigate the optimum cavity depth. 

For each porosity, the static pressure distributions, Schlleren 
photographs and Impact pressure wake surveys were obtained over a range 
of flow Mach number of 0.65 to 0.82. 
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1. Supercritical airfoil 

2. Porous surface 

3. Cavity beneath porous surface 

4. Free stream conditions 

5. Embedded supersonic region 

6. Sonic line 

7. Terminating shock wave 

8. Flow circulation through the porous surface 

9. Wake survey rake 
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Fig. 9 Passive drag control for supercritical airfoil 
at transonic Mach numbers 
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Fig. 10 Conventional and supercritical airfoils 
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CALIBRATION OF THE EXPERIMENTAL FACILITY 

For Che Initial calibration of the Transonic Wind Tunnel with 
Che constant area test section, as shown in Figs. 1 and 3a, axial static 
pressure distributions and Schlleren and shadoi^raph photographs were 
obtained. Different airfoil models were investigated in Che original 
constant area mahogany test section, and from these preliminary test 
results it was decided Co open Che Cop, bottom, and side walls as shown 
in Fig. 3b. Experiments were conducted in the modified test section with 
the double wedge and biconvex airfoils mounted in the center of the test 
section, as well as with the 3-inch chord circular arc airfoil mounted 
on the bottom surface as shown in Fig. 6. The results for both Che orig- 
inal and the modified test sections with and without the airfoils in the 
test section are discussed below. 

3.1 Constant Area Test Section 

3.1.1 Static Pressure Variation with Test Time 

The variation of the geometrical area ratio A/A*, with the 
axial distance is presented in Fig. 12, where A* is the area of the orig- 
inal constant area test section and A is Che area corresponding to the 
axial distance of the wind tunnel. A large area ratio for the settling 
section was selected to keep a low velocity through the entrance, honey- 
comb, and screens with corresponding low pressure loss. The large con- 
traction ratio before the test section decreased the turbulence level 
in the accelerating flow region. To recover as much total pressure as 
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possible downstream of the test section, the area ratio at the end of 
Che diffuser section was equal to 5, which Is the maximum available with 
the existing 16-lnch diameter discharge piping system. 

To Investigate the steadiness of the flow field and the per- 
formance of the different components of the tunnel as well as the test 
time, the static pressures at selected locations throughout the tunnel 
and the vacuum tank pressure were taken for various flow durations. The 
time variations of the static pressure at different axial locations are 
presented In Fig. 13. The static pressures In the settling chamber, 

X * 24 Inches, and at the end of the test section, x > 109.25 inches, 
were nearly constant with time. Indicating steady flow, but In the dif- 
fuser, X * 122 Inches, downstream of the test section, the static pres- 
sure increased with time as the pressure in the large vacuum tank built 
up. 

With the initial test section area, the steady flow duration 

in the test section was approximately 9 seconds. This was confirmed 

later with the aluminum test section by observing the steadiness of the 

shock wave through a frosted glass in the back of the camera box. A 

20 

Klstler quartz pressure transducer with a response time of 10 micro- 
seconds was Installed In the test section to obtain the Instantaneous 
variation of the static pressure with time. The pressure gage output 
was recorded on an oscilloscope and the result indicated that the flow 
In the test section was established In 0.4 second and remained constant 
approximately 9 seconds. This test time was long compared to the test 
time of 8 milliseconds for the Ludwleg tube used to obtain the pressure 
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dlscrlbutions and the Interferogram photographs at Mach numbers of 0.7 
and 0.8 over a 10-percent biconvex airfoils with air, SF^ and CO 2 gases 
in Ref. 21. 

3.1.2 Axial Static Pressure and Mach Number Distributions 

Based upon the nozzle cross sectional area distribution for 
the original test section presented in Fig. 12, the corresponding theo- 
retical static pressure ratio and Mach number distributions are pre- 
sented in Figs. 14 and 15, respectively. In the calculations of these 
flow parameters, the flow was assumed to be adiabatic, Invlscid, one 
dimensional, and steady as discussed in Ref. 22. Subsonic and super- 
sonic flow conditions were calculated downstream of the test section. 

The nozzle static pressure ratio distributions for different 
flow duration time are compared with the theory in Fig. 14. For a flow 
duration of 2 seconds, the supersonic flow remained attached to the 
diffuser walls at least to a distance of 11 inches downstream of the 
test section. This was possible because the vacuum tank pressure for 
this flow duration was very low, as shown in Fig. 13. As the flow dura- 
tion Increased and brought about a corresponding increase in the vacuum 
tank pressure, the region of the supersonic flow downstream of the test 
section became shorter. At a flow duration of 8 seconds, the flow became 
subsonic slightly downstream of the test section, and the static pressure 
in the diffuser approached the theoretical subsonic flow values. 

The static pressure distribution from the settling chamber to 
the test section exit remained nearly the same for flow duration of 4 to 
8 seconds. The static pressures for the flow duration of 2 seconds were 
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slightly higher than the values for longer test times. This difference 
was caused mainly by the long response time of the mercury manometers 
with long plastic tubings. A quartz piezoelectric pressure transducer 
was used for measuring the pressure In the test section to determine the 
flow establishment time as discussed previously. 

The static pressure distributions in the test section were 
much higher than the theoretical Invlscld values. Because of the boundary 
layer growth along the test section walls with the original constant area, 
Fig. 3a, the flow was choked near the exit of the test section i^ere the 
effective flow area was a minimum. Ahead of the choked sonic Mach number 
location, the flow was subsonic and the static pressure ratio was greater 
than the sonic value of 0.528. The test section was modified to correct 
for this boundary layer choking phenomenon. Fig. 3b, and the results will 
be discussed later. 

The axial Mach number distributions for different flow durations 
are compared with the Invlscld theory In Fig. 15. For a test time of 2 
seconds with low vacuum tank pressure, the flow accelerated to a Mach 
number of approximately 1.72 in the diffuser before the flow separated 
from the walls. Again as the flow duration Increased, the maximum super- 
sonic flow velocity and the region of the supersonic flow In the diffuser 
downstream of the test section decreased. The Mach number dlstlrubtions 
from the settling chamber to the sonic location In the test section were 
nearly Identical for flow durations of 4 to 8 seconds. This axial Mach 
number variation was due to the boundary layer growth along the tunnel 
walls as discussed previously. 
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3.2 Modified Test Section 

The axial static pressure distributions* Fig. 14, for the 
original constant area test section, Indicated the choking of the flow 
near the exit of the test section due to the boundary layer growth along 
the tunnel walls. To correct for the boundary layer growth, the follow- 
Ing test section modifications, as discussed earlier, were Investigated: 
opened top wall; opened top and bottom walls; opened top, bottom, and 
side walls, as shown In Figs. 3a and 3b. The Mach number distributions 
In the test section for these modifications are presented In Fig. 16. 

The modifications of the test section and the results are discussed 
below. 


3.2.1 Constant Area Test Section 

The axial distance at the entrance of the test section was 
73.5 Inches, i^lle the exit was located at 110.0 Inches. With the orig- 
inal 3- Inch by 15.4-lnch constant area test section, the Mach number In- 
creased from 0.72 at the entrance to 1.10 at the exit due to the boundary 
layer growth along the walls, as shown In Fig. 16. The flow was sonic a 
few Inches before the test section exit and became supersonic at the exit 
due to the presence of the diverging diffuser walls, which tended to thin 
the boundary layer at the test section exit. 

3.2.2 Test Section with Opened Top Wall 

In order to accommodate for the boundary layer growth, the top 
wall of the test section was opened by tapering the mahogany wall from 
15.4 Inch opening at the entrance to 16.1 Inches at the exit. This modi- 
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flcaclon made tha Mach number dlacrlbutlon In the teat aectlon flatter 
than that of the original configuration as shown in Fig. 16. The Mach 
number varied from 0.75 at the entrance to 1.02 at the exit. In addition, 
the sonic location was closer to the exit than In the original constant 
area test section. 

3.2.3 Test Section with Opened Top and Bottom Walls 

The bottom wall was tapered over the rear 16.5 Inches from 
15.4 Inches to 16.1 Inches as shown In Fig. 3b. The Initial 20 Inches 
of the bottom wall was horizontal to permit the investigation of the 
shock wave/boundary layer interaction over circular arc and supercritical 
airfoils with a porous surface. The Mach number distribution with the 
top and bottom walls opened was very flat over the first 20 Inches of 
the test section as shown in Fig. 16. Along the test section, the Mach 
number Increased from 0.88 to 0.98 at the exit. For this configuration, 
the flow was choked near the diffuser entrance. 

3.2.4 Test Section with Opened Top. Bottom and Side Walls 

For further Increase of the free stream Mach number, the test 
section area was enlarged by tapering the side walls from a width of 3.0 
Inches at 20 Inches downstream of the entrance to the test section to a 
width of 3.30 inches at the exit of the test section as shown in Fig. 3b. 
The Mach number increased from 0.89 at the entrance to 1.26 at 30 Inches 
downstream and then decreased to 1.13 at the exit. The flow was sonic 
at approximately 17.5 inches from the entrance. By opening the 4 walls 
of the test section, the Mach number at the entrance was Increased from 
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0.72 for the original constant area test section to 0.89. The use of 
mahogany for the top and bottom walls and thick Incite for the side walls 
had made the modification of the test section to correct for the boundary 
layer growth fairly easy. Based on the knowledge gained from the cali- 
bration of this Initial mahogany test section, an aluminum porous top wall 
test section was designed and constructed. This test section helped In 
obtaining the flow data without side wall deflection effects, besides 
being more rigid and durable. 

3.3 Boundary Layer Surveys for Test Section Walls 

The uniformity of the flow field and the choking Mach number 
in the test section depended on the growth of the boundary layer along 
the walls. By correcting for the boundary layer growth. It was possible 
to control both of these flow phenomena. For this purpose, boundary 
layer surveys were made at selected locations on the bottom and side 
walls of the original constant area test section. Fig. 3a. Since the 
test section was symmetrical, the boundary layer growth was assumed to 
be the same on the top and bottom surfaces and also on both side walls. 
Boundary layer displacement thicknesses were calculated for both walls 
from the liiq>act pressure surveys across the boundary layers. 

3.3.1 Bottom Wall Boundary Laver Survey 

Locations near the Inlet, middle and exit of the test section 
at corresponding axial distances of x ■ 76.75, 94.87 and 106.75 Inches 
were selected for the boundary layer survey for the bottom wall. Impact 
pressures across the boundary layer were measured with a small Impact 
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pressure probe and the results are presented In Fig. 17a. Inpact pressure 
surveys Indicated that the boundary vas fairly thin at the entrance to 
the constant area test section and becane quite thick at the exit. 

From the measured static and Impact pressure data for the bound- 
ary layer at different axial locations In the test section* the correspon- 
ding velocity distributions In the boundary layer were calculated for an 
assumed adiabatic flow. The velocity In the boundary layer was normalized 
with respect to the free stream velocity. The vertical distance from the 
wall was normalized with respect to the boundary layer thickness. The 
normalized velocity distributions for the three axial locations are pre- 
sented In Fig. 17b. At all three locations the velocity distribution 
Indicated that a turbulent boundary layer existed over the surface. 

3.3.2 Side Wall Boundary Layer Survey 

The side wall boundary layer was surveyed on the centerline of 
the test section at three axial locations of 77* 92.62 and 106.50 Inches* 
idiich corresponded to near the Inlet* middle and near the exit of the 
test section* respectively. The results of these Impact pressure surveys 
are presented In Fig. 18a. The side wall boundary layer Increased with 
the axial distance, and at the test section exit* x ■ 106.5 Inches* the | 

boundary layer was fairly thick. To correct for this boundary growth* j 

the side walls were widened as Indicated In Fig. 3b. I 

The normalized velocity distributions calculated from the mea- | 

sured static and Impact pressures for the three axial location In the | 

1 

test section are presented In Fig. 18b. The side wall boundary layer 
distributions for axial locations of 77 and 106.5 Inches are ccjmpared | 

i 
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wlch ch« laminar and turbulent boundary layer velocity dlatrlbutlons in 
Fig. 18c, and the reaulta Indicate the turbulent boundary layer char- 
acterlatlca. 


3.3.3 Boundary Laver and Dlanlacement Thlckneea on Teat Section Welle 
The boundary layer thlckneea for the teat aectlon bottom and 
aide walla, determined from the impact preeaure aurveya at three axial 
locatlona, Flga. 17a and 18a, are preaented In Fig. 19a. It ia evident 
In thla figure that the boundary layer thlclmeae and the growth along 
the bottom and aide walla were about the aaa». The apparent boundary 
layer thlckneea near the exit on the aide wall waa allghtly leaa than on 
the bottom wall becauae the liapact preat «re aurvey waa not extended 
further from the aide waU., aa ahown in Fig. 18a. The boundary layera 
on the tunnel walla grew more rapidly over the front part of the conatant 
area teat aectlon. Fig. 3a, than over the aft portion. 

From the meaaured Impact preaaure acroaa the boundary layer 
and the wall static preaaure, the velocity and denalty were calculated 
by aaaumlng an adiabatic flow. Theae valuea were uaed to calculate the 
boundary layer dlaplacement, 6*, from the equation 
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where u ia the velocity, p the denalty, and p^u^ ia the mesa flow rate 
outside the boundary layer. The calculated dlaplacement thlckneea for 
the bottom and aide walla are preaented in Fig. 19b. For both walla. 


the displacement thlckneea grew rapidly over the Initial portion of the 
constant area test aectlon. The displacement thickness on the bottom 













wall near the exit of the test section was less than the value at the 
axial location of 94.87 Inches. This decrease In the displacement 
thickness was caused by the acceleration of the flow toward the exit of 
the constant area test section as shown In Fig. 16. 

3.4 Flow over Double Wedge. Biconvex and Circular Arc Airfoils 

3.4.1 Flow over Double Wedge Airfoil 

A 10-percent-thlck double wedge airfoil with a chord of 3 
Inches was constructed out of aluminum to Investigate the flow field and 
the shock wave structure at transonic Mach numbers. The airfoil was 
mounted In the center of the plexiglass side walls test section at an 
axial distance of 91.75 inches. Static pressure orifices are located 
on the top wedge surface at 40, SO and 60 percent of the chord from the 
model leading edge. 

For these tests, the top, bottom and side walls of the tunnel 
were opened as shown In Fig. 3b. From the measured static pressures and 
the reservoir total pressure i the local flow Mach number in the test sec- 
tion and over the model was determined. Ahead of the airfoil the flow 
Mach number was 0.88. Over the wedge the flow accelerated to a Mach 
number of approximately 1.26 at a location slightly downstream of the 
mid-chord as shown in Fig. 20a. Only a few static pressure orifices were 
Installed In this airfoil to obtain the maximum flow velocity over the 
wedge surface with a simple model. 

The Schlleren and shadowgraph photographs of the flow over the 
10-percent thick double wedge airfoil were taken to observe the shock 



waves and boundary layer in the transonic flow region over the airfoil. 
The Schlleren photographs with vertical and horizontal knife edge 
positions are presented in Figs. 20b and c. In the vertical position, 
the knife edge was more parallel to the shock waves and made the shock 
wave and the expansion region more visible, as shown in Fig. 20b. The 
expansion fan from the maximum thickness region of the airfoil is visi- 
ble as a dark or light region depending on how the vertical knife edge 
cut the light beam from left or right. In this region, the density de- 
creased as the flow accelerated through the fan. An oblique shock wave 
terminated the expansion region and Intersected the nearly normal ter- 
minating shock wave. In the photographs the stress concentration in the 
lucite side walls around the model support holes is visible as well as 
the plastic tubing outside the walls, used for measuring the static 
pressures on the model. 

The Schlleren photograph, with the horizontal knife edge to 
make the boundary layer over the airfoil and the wake more visible, is 
presented in Fig. 20c. In this photograph, the boundary layer near the 
trailing edge and the wake flow are quite visible, but only the oblique 
shock wave is distinctly visible. 

3.4.2 Flow over Biconvex Airfoil 

An aluminum, 12-percent- thick, biconvex airfoil with a 3-lnch 
chord was placed in the center of the test section, similar to the double 
wedge airfoil. Static pressure orifices on the top of the airfoil were 
located at x/c of 0.21, 0.39, 0.50, 0.62 and 0.79. The Mach number ahead 
of the model was 0.89. The flow accelerated over the biconvex airfoil to 
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Fig, 20b Schlieren photograph of flow over 10-percent-thick double 
wedge airfoil, vertical knife edge ,M*0,83 


Fig. 20c Schlieren photograph of flow over 10-percent-thick double 
wedge airfoil , horizontal knife edge, M-0.83 







a maximum valxte of 1.25 downscream of the mid-chord as shown In Fig. 21a. 
The Schlleren photograph of the flow over the 12-percent-thlck biconvex 
airfoil shows a nearly normal shock wave, which is more visible on the 
upper surface. The shock wave location on the upper surface was at 78 
percent of the chord, while on the lower surface. It was slightly down- 
stream of this location. The photograph shows also some distortion of 
the airfoil leading edge, because of the side walls deflection. This 
shortcoming was overcome later. In the aluminum test section by going to 
a thicker side wall plexiglass and adding aluminum plates on both sides 
with a circular cut at the model location to be able to take the flow 
field and the model photographs. From the shadowgraph photograph. Fig. 
21b, the boundary layer and Its separation over the 12-percent-thlck 
biconvex airfoil are visible. 

3.4.3 Flow over Circular Arc Airfoil 

A 12-percent-thlck convex airfoil made out of plexiglass was 
mounted on the test section bottom surface, where the concept of the 
shock wave/boundary layer Interaction control over a 14-percent- thick 
supercritical airfoil was Investigated. Six static pressure orifices 
were used to obtain the pressure distribution over the model surface, 
and the Mach number distribution Is presented In Fig. 22a. The Schlleren 
photograph showing the shock wave behavior Is presented In Fig. 22b. 

With the vertical knife edge, the termination shock wave Is very distinct 
In the photograph, due to a maximum flow Mach number of 1.29 over the 
model. At this high transonic Mach number, the pressure and the density 
change across the shock wave was appreciable, as indicated In Ref. 22. 
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The relatively strong shock wave Indicated the presence of a large super- 

23 

sonic flow region over the model. 
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Fig. 21b 


Shadowgraph photograph of flow over 
-thick biconvex airf oil ,.M»0. 83 
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PART 4 
THEORY 

The basic equations of compressible flow used In the present 
research program for the data reduction are as presented In this section 
for the assumption of steady, perfect gas flow. 

4.1 Isentroplc Flow Equations 
Determination of Mach number 

The Mach number In the teat section, over the model surface 
and throughout the tunnel, was obtained from measurements of the static 
pressure, P, and the total pressure In the settling chaad>er, P^. The 
Mach number Is then given by the well known compressible Isentroplc flow 
relations 

. {1 + illil Y-1 ^ (1) 

^ere y ■ taken as 1.400. 

Determination of Temperature 

By reading the total temperature, T^, In the settling chamber, 
which turns out to be very close to the room temperature, the temperature, 
T, at any location In the tunnel can be determined, once the corresponding 
Mach number Is known, through the adiabatic relation: 

- (1 + ^ M^)“^ (2) 

o 

Determination of Density 

Using the perfect gas relation to get the total density p^. 
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The density et any tunnel locationy for a corresponding Mach nuaber, can 
be derived from the Isentroplc, perfect gas relation: 

_ _ 1 _ 

^ ■ (1 + M^) (4) 

^o 

where P and T are the total pressure and total teo^erature measured In 
0 0 

the settling chamber. 

Determination of Speed of Sound 

Knowing the Mach number » the speed of sound can be derived 

from the adiabatic relation and is given by 

1 

" ”2 

- (1 + ^ M^) , (5) 

idiere a^ Is the speed of sound of the gas at rest. By measuring in 
the settling chamber, a^ Is derived from the speed of sound relation 

a^ ■ /yRIq , (6) 

and the local velocity sound, a, can also be determined by the same 
relation 

a - /yH (7) 

once the corresponding T Is derived from Eq. (2). 

This speed of sound relation can be approximated by 

a * 49»^ ft/sec (8) 

for air If T Is In degrees Ranklne. 
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Determination of the Flow Velocity 

Knowing the Hach number and the speed of sound, the corre- 
sponding local velocity can be determined from the Mach number relation- 
ship 


M - 7 (9) 

a 

Determloatlon oi: the Pressure Coefficient 

The pressure coefficient In compressible flow depends on the 
Mach number as well as the static pressure and is defined as 


C 

P 


P-P 


1 TT 2 


( 10 ) 


For an Isentropic flow, the equation may be rewritten in terms of the 
local Mach number as 


2 + (y-1)M ‘ Y-1 


[- 


ym„ 


2 + (y-I)M^ 


] 


- 1 


( 11 ) 


M Is the free stream Mach nximber measured ahead of the model and M Is 

00 

the local Mach number obtained from the Eq. (1) by measuring the local 
static pressure and the total pressure in the settling chamber. The 
steps between Eqs. (10) and (11) can be found In Ref. 24. 

The Area-41ach Number Relation 

The area-Mach number relation for a steady, Invlscld, compressi- 
ble and one-dimensional flow is given by 


2 ^ 
h 


( 12 ) 


where A* Is the throat area. 
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4.2 Normal and Oblique Shock Wave Relations 
Normal Shock Wave Relations 

The relation between the Mach numbers upstream and downstream 
of the normal shock wave Is given by 



1 + 


I2i 

2 





(13) 


where is the Mach number downstream of the shock wave and M^ is the 
Mach number upstream. Other useful relations can be derived from the 
continuity equation, as the ratio of densities 

p U. (y+1)M,^ 

7^ “ TT T — 

•^1 ^2 ( y - I ) Mj ^ +2 

The ratio of the static pressures is often used to define the shock wave 
strength and is given by 


From the energy equation the temperature jump across the shock wave is 
given by 


i . llrtil 2 . 1, (16) 

^1 (Y+1) M^^ ^ 

The corresponding jumps in density, pressure and temperature 
are from lower to higher values across the shock wave. The shock is said 
to compress the flow. 

An Important result is the increase in the entropy through 


the shock wave given by 
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(i.7) 


frtilch means chat Che increase in entropy is third order in shock strength. 
The change in entropy can also be related to the total pressures upstream 
and ''.j«mstream o£ the shock wave: 


"2-^1 


R 


In 


^ol 


(18) 


^o2 


fdiere the ratio of total pressures is given by 


Pol 2 yMj^^-(Y-1) 


^ Y+DM 2 


(19) 


( y - 1 ) Mj ^‘+2 


Oblique Shock Wave Relations 

The same relations used for the normal shock wave can be 


applied here for the oblique shock wave with the modification of and 


being replaced by their normal components Musing and M 2 sin( 6 - 6 ), re- 
spectively. B being the flow inclination to the shock wave and 6 is the 
flow deflection after the oblique shock wave, 



. . Shock wave 


This gives the corresponding relations for the oblique shock: 

1 


M2^sin^(e-e) 


^ Iri w 2 . 2„ 
+ -‘rr- M, sin B 


YMj^^sin^e - 


( 20 ) 


(Y+l)Mj^^sin^6 


(Y-DMj^^sin^e+2 


( 21 ) 






^2 2(y- 1) “l 2 2 

^ . 1 + liJLJi _L^ ^(yM. ^ 8in^0+l) 

^1 (Y+1) 6 


S 2 -S 1 


r^- l- 


^ol 2 yM, ^ sin^g-(Y-l) (y- 1)M, ^sin^g+2 


(Y+l)M.^8in^e ^ 
\ 22 


4.3 Profile Drag Derivation 


The airfoil aection profile drag meaaurementa were computed 
from the wake survey rake measurements by the method of Ref. 25 utilizing 


the following equations: 
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d c 


ri 1 

Cl - 2(^) i' (f-)^ 


I Y- 


Izi 1 

Y \ 2 


P izi 

l-(f ) ^ 

00 


l-(^) 


l-(^) 

00 


1=1 1 

^ Y \2 


with » free stream total pressure 

a free stream static pressure 

U > local total pressure in the wake 
P » local static pressure in the wake 
s point drag coefficient 

C, = section drag coefficient 
d 

To obtain the section drag coefficients, point drag coefficients were 


computed for each total pressure measurement in the wake by using Eq. (26). 
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These point drag coefficients were then summed up according to Eq. (25) » 
by numerical integration across the wake, based on the trapezoidal 
method. 


4.4 Boundary Layer Analysis 

Lam<nar Boundary Layer Along a Flat Plate 

Starting from the boundary layer equations for a steady state 
and incompressible flow and introducing new dimensionless variables we 
get the following Blaslus equation: 


f f»* + 2f" - 0 

where f is a dimensionless stream function defined as: 


dn u 


(27) 


(28) 


with 
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(29) 


a numerical solution of Eq. (27) by L. Howarth is presented in Ref. 26. 

The boundary layer thickness deduced from the numerical solu- 
tion is given by 

. 0 /^ 


6 » 5 

and based on the Reynolds number, R 


(30) 
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■ , and the boundary 


layer thickness can be written as 


5. Ox 
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( 31 ) 


the dimensionless coordinate n can be written as 
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vx X e„ 


which yields 


6 5.0 


(32) 


and the values of y/6 are presented in Ref. 26. 

A physically meaningful measure for the boundary layer thick- 
ness is the displacement thickness, The displacement thickness is 
that distance by which the external potential flow field Is displaced 
outwards as a consequence of the decrease In velocity In the boundary 
layer. The decrease In volume flow due to the Influence of friction Is 
given by 

/ (u^-u)dy 
'y»o 

so that for 6 * we have the definition 
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From the solution presented In Ref. 26, the displacement thickness 6* 
is given by 
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Turbulent Boundary Layer Along a Flat Plate 


The velocity profile of a turbulent boundary layer for an In- 


compressible flow can be represented by an empirical equation' 

1 
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(36) 


and the assumption of a simple -^h-power law agrees with experiment. 


Based on -yth-power law which Is valid for moderate Reynolds numbers, the 
velocity profile of a turbulent boundary layer over a flat plate at zero 
Incidence is given by 


^ - (4)" 


(36a) 


where 6 « 5(x) denotes the boundary layer thickness which Is a function 
of distance, x. The velocity distribution Is presented In Figs. 18c and 
19c. The equation for the local shearing stress Is given by 
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where 6 is the momenr*im thickness. 

Using equation (36a), and the definition of 6 In Eq. (37), we 

obtain 


6 * - g 
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( 38 ) 
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The differential equation for the boundary layer thickness Is given by 
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The integration of Eq. (39) yields 
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*x I*. 

The displacement thickness, 6*, from Eq. (38) is given by 

... _ .046x 


Compressibility Effect 


The above equations have been derived for an incompressible 
flow case. For a laminar boundary layer, Fig. 15.6 in Ref. 26 shows that 
for a free stream Mach number ranging from 0 to 1.00, which corresponds 
to the present experimental work, the compressibility effect? on the 
velocity distribution within the boundary layer over a flat plate of 
zero Incidence are negligible. Therefore, a comparison of the experi- 
mental boundary layer data with the theory of incompressible flow is a 
good approximation for our objective. 
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DISCUSSION OF THE CALIBRATION PHASE RESULTS 
5.1 Transonic Wind Tunnel 

A 3-inch X 15.4-lnch transonic wind tunnel was constructed 
(as shown in Figs. 1 and 2} which utilized the existing dryer, vacuum 
tank and vacuum pumps for the 4-Inch by 4-Inch Mach 3 Supersonic Wind 
Tunnel. The calibration of the flow through the tunnel was determined 
from the static pressure distribution throughout the length of the tunnel, 
the Impact pressure in the settling chamber, and the impact pressure sur- 
veys in the test section. Schlleren and shadowgraph photographs were 
taken for a lO-percent-thlck double wedge and 12-percent-thlck biconvex 
airfoil mounted in the center of the test section, as well as for a 12- 
percent-thick circular arc airfoil placed on the test section bottom 
surface. The flow measurements were taken for different flow durations 
to determine the time of the flow establishment and the duration of the 
transonic flow test time. 

The honeycomb and screens, Fig. 1, being the best combination 
for reducing turbulence at the entrance to the settling chamber, as dis- 
cussed in Ref. 27, caused the flow entering the tunnel to be steady, as 
indicated by the constant settling chamber impact pressure shown in Fig. 
13. The multilayered aircraft pl 3 rwood with external reinforcement used 
for the contraction section, see Figs. 1 and 2, performed well for all 
test section Mach numbers. A grad\ial contraction profile accelerates 
the flow which becomes uniform at the entrance to the test section. The 
thick luclte side walls for the test section made it possible to observe 
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Che shock waves and boundary layer over the airfoil models mounted In 
the center and on the bottom wall of the test section. 

The long gradual diffuser section downstream of the test sec- 
tion performed well as indicated by Che good pressure recovery in the 
diffuser, shown in Fig. 14. This permitted a choked, sonic velocity at 
Che end of the original constant area test section. Fig. 3a, %d.th a small 
pressure drop across the transonic tunnel, and allowed for a longer test 
time with Che available vacuum tank and pumps. 

The pneumatic, 16-lnch diameter, quick-acting valve located at 
the exit of the diffuser section. Figs. 1 and 2, opened and closed in 
approximately 1.5 seconds with a good vacuum seal. The piezoelectric 
pressure measurement in the test section Indicated that the flow was es- 
tablished In 0.4 second. 

5.2 Instrumentation 

The technique for Che pressure data acquisition with the mano- 
meter board. Fig. 2a, was Improved by replacing the ball valves with an 
electrically timed and actuated solenoid valve system. Leakage tests 
of up to 8 hours In length were conducted with negligible pressure losses 
In Che measuring system. The system test results Indicated that the time 
response was sufficiently short and accurate to obtain reliable pressure 
data, as did the fact that all of the pressures were taken simultaneously, 
which was not the case with Che ball valves. 

From Figs. 20b and 22b, for flows over a 10-percent-thick double 
wedge, and 12-percent- thick convex airfoils, respectively, it was apparent 
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that the Schlleren system was capable of obtaining photographs of the 
shock waves and boundary layers over these models. To observe the shock 
waves and boundary layer, the knife edge was positioned In the vertical 
and horizontal positions, respectively. 

As can be seen from the shadowgraph photograph shown In Fig. 

21b, the shadowgraph system appeared to be functioning satisfactorily. 

Further Improvements of the optical system were made by the use of thicker 

luclte side walls, reinforced by aluminum plates, to eliminate the wall 

13 

deflection on the new aluminum test section. 

5.3 Tunnel Wall Boundary Laver 

As expected, the thickness of the boundary layer over the bottom 
and side walls for the constant area test section Fig. 3a, was Increasing 
as It progressed from the entrance of the test section to the exit, as 
seen In Fig. 19a. For all axial locations, the velocity profiles for 
both walls agreed closely with the fully developed, turbulent boundary 
layer profile, as shown In Figs. 18c and 19c. On the bottom wall the 
boundary layer thickness Increased from 0.094 Inch to 0.53 Inch near the 
exit. This relatively thick boundary layer over the test section walls 
was the reason for modifying the original constant area test section, as 
shown In Figs. 3a and 3b. Data from the liupact pressure surveys across 
the boundary layer, were used to calculate the boundary layer displace- 
ment thickness for the side and bottom walls. The results are presented 
In Fig. 19b. The displacement thickness Increased steeply over the 
Initial region of the test section, but Its growth from the center to 
the exit was more gradual. The displacement thickness was slightly 
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greater on the bottom wall, but near the test section exit, the thick- 
nesses were about the same. 

5.4 Flow Conditions in the Original and Modified Test Sections 

With the original constant area test section shown In Figs. 1 
and 3a, the Mach number at the entrance was 0.714 und Increased monotoul- 
cally over the 36.5 inch length to a Mach number of 1.11, as shown In 
Fig. 16. This Increase In the Mach number In the test section was caused 
primarily by the growth of the boundary layer along the tunnel walls. The 
sonic velocity occurred close to the exit of the test section. To correct 
for the boundary layer displacement thickness, the top wall of the test 
section was opened 0.7 inch at the exit. Fig. 3b. The entrance Mach number 
Increased to 0.75 and the exit Mach number decreased to 1.08, as shown In 
Fig. 16. By opening the top and bottom walls. Fig. 3b, the Mach number 
attained a more uniform value of 0.88 throughout most of the test section, 
except at the exit where the Mach number Increased to 0.98, as shown In 
Fig. 16. 

A third test section modification was made to Increase the Mach 
number at the entrance to the test area. This was accomplished by keeping 
the side walls parallel over the initial 20 Inches and then tapering them 
to 0.15 inch at the exit. Fig. 3b. The Mach number at the entrance In- 
creased to 0.89 and reached scnlc velocity slightly ahead of the beginning 
of the diverging side walls. Fig. 16. The flow accelerated to a peak Mach 
number of 1.26 in the diverging side wall section, and subsequently de- 
creased to a Mach number of 1.14 at the exit. The addition of the side 
wall modification appeared to have made the geometric area Increase faster 
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Chan Che boundary layer grouch along Che cesc secclon walla. Thus, che 
flow acceleraced co a peak Mach number of 1.26 before che Increase in 
boundary layer grouch became greacer chan che Increase In Che geomecrlc 
area. 

Following che resulcs of chese callbraclon ceacs, an aluminum 

cesc seccion wich an adj useable cop porous wall was designed and con- 
13 

scrucced. This new Cesc secclon used many of che f secures from ics 

mahogany predecessor, buc ac che same clme, modlflcaCions were made co 

overcome some previously observed deficiencies. The luclce side walls 

were made chlcker and aluminum places were added Co che sides co mlninilae 

che wall deflecclon. The variable porous cop wall helped co reduce Co a 

minimum che wall interference and Che model blockage efface ac high sub- 

28 

sonic Mach numbers. Ocher mechods of minimizing che wall InCerference 
effaces use slocced and adapclve walls as discussed in Refs. 29, 30 and 31. 

5.5 Flow over Double Wedge and Biconvex Airfoils 

The 10-percent- Chick double wedge airfoil wiCh a 3-inch chord 
was mounced in Che cenCer of the test section. With a Mach number of 
0.88 ahead of che model, the maximum local Mach number over Che airfoil 
reached 1.26. The Schlleren photograph. Fig. 20b, indicates the flow ex- 
pansion region and the oblique and termio&tlon shock waves. The nearly 
normal cermlnaclon shock wave was located at 68 percent of the chord. 

The Schlleren photograph and Che pressure data are similar co Che results 
obtained in Ref. 14 for a 10-percent-thick double wedge 3-inch chord air- 
foil in a 3-lnch by 16- inch transonic wind tunnel wich solid walls. 
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For the 12-perceat-chlck biconvex airfoil with a 3-lnch chord 
and a free stream Mach number of 0.89, the flow accelerated to a peak 
Hach number of 1.25 downstream of the mid-chord, before decelerating 
towards the trailing edge. The Schlieren photograph shows a nearly 
normal shock wave at 78 percent of the chord on the upper surface of 
the model, while on the lower surface it is slightly downstream. This 
difference in locations comes from the fact that the angle of attack was 
not quite tero. On the shadowgraph photograph. Fig. 21b, the boundary 
layer and its separation are visible. In Ref. 15, the Hach nmnber dis- 
tributions over a 10-percent- thick biconvex airfoil were obtained over 
a Mach number range of 0.80 to 0.90 in a solid test section. In this 

reference the actual Mach number distribution at M of 0.89 is seen to 

00 

have a higher local maximum Mach number. This is because of the differ- 
ence in the thickness between the two models. In Ref. 16, strong shock 
waves are observed over a 12-percent-thlck biconvex 3-lnch chord air- 
foil, for a free stream Mach number of 0.89. The location of these shock 
waves is the same as the one obtained over the present similar biconvex 
airfoil. The tests in Ref. 16 were conducted in a 2-inch by 20-inch, 
solid wall, transonic wind tunnel. 

5.6 Flow over Circular Arc Airfoil 

The 12-percent- thick circular arc airfoil with a 3-lnch chord 
was placed on the bottom wall of the test section. For a free stream 
Mach number of 0.88, the maximum local Mach number obtained over the model 
was 1.29. The model had a similar Mach number distribution as the 10- 
percent-thlck model in Ref. 15, with a slightly higher maximum local 


Mach number, because of Che difference in Che model Chickneas. The large 
excenslon of the scrong shock wave Inco Che flow Indlcaces Che presence 
of a large supersonic zone over che airfoil surface, as shown In Che 
Schlleren photograph of Fig. 22b. 

5.7 Passive Shock Wave-Boundary Laver Control for Supercritical Airfoils 

To invesclgate Che concepc of che passive shock wave-boundary 

layer concrol for possibly reducing the drag of supercritical airfoils, 

various surface porosities were tested with different cavity geometries. 

The preliminary tests of the effects of porosity and cavity size on che 

shock wave/boundary layer interaction were conducted on a 12-percent-Chlck 

circular arc airfoil with a 4-lnch chord to determine the effective porous 

surface configuration. Based on these test results, a 14-percent- thick 

18 

NASA supercritical airfoil was investigated with selected surface poro- 
sities and cavity sizes. 
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PART 6 


EXPERIMENTAL RESULTS AND DISCUSSION 

The following experiments have been conducted In the aluminum 
porous top wall test section described previously. A 12-percent- thick 
circular arc profile and a l4-percent- thick supercritical airfoil, both 
with a 4-inch chord, were investigated on the bottom surface of the test 
section. Each model was mounted in the middle of the first 20-lnch flat 
surface. Two adjustable wedges Installed downstream of the model at the 
test section exit. Fig. 4, were used to control the free stream Mach 
number from low subsonic to transonic speeds. The major part of the ex- 
periment was conducted with no test section top wall porosity, and a few 
tests were made to study the wall Interference on the flow field. 

The experiments which will be described below are aimed at re- 
ducing the drag by controlling the shock wave/boundary layer Interaction 
over conventional and supercritical airfoils. The circular arc and super- 
critical airfoils were selected to achieve the desired objective. Both 
airfoils were Investigated over a range of subsonic and transonic Mach 
numbers with various model porosities and different cavity sizes located 
in the region of the shock wave/boundary layer Interaction. 

The way to achieve passive drag control, as discussed previously. 
Is to install a porous surface with a cavity underneath at the location 
of the shock wave over the model. By this procedure, because of the pres- 
sure difference across the shock wave, a part of the boundary layer flow 
will be moved from the region doimstream to the region upstream of the 
shock wave through the porous surface and the cavity, as shown In Fig. 9. 
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This method will Chin Che boundary layer, downstream of the shock wave, 
by removing the decelerated flow, thus minimizing Che separation. It 
also produces compression waves ahead of the shock ^Ich decreases Its 
strength. The combination of these phenomena helps In reducing the drag 
of Che airfoils In transonic flows, which Is the major goal of the current 
aerodynamic research. Therefore, the research described In this Section 
Is aimed at Increasing the technology base for the development of practi- 
cal means of Improving aerodynamic performance and efficiency at high sub- 
sonic speeds. 

Before Initiating the airfoils Investigation, Che aluminum test 
section was calibrated with and without Che side wall aluminum plates 
added to the luclte side walls. For the same control wedge setting, the 
test section with the aluminum plates gives a higher free stream Mach 
number than without because the deflection of the luclte side walls In 
the absence of the aluminum plates makes Che effective min limim area smaller 
near the middle of the test section. By assuming the free scream area 
ahead of the model Co be nearly constant and less affected by the deflec- 
tion of the side walls, the ratio of the free stream area over the mini- 
mum area Is smaller with the aluminum plates. This leads to a higher 
free scream Mach number for a subsonic speed range, as shown by the area- 
Mach number relation, Eq. 12, presented In Fig. 23 and confirmed by the 
experiments, Fi^re 24 • Consequently, the Mach number dis- 
tribution Is higher. The test results also show a more uniform free 
stream Mach number distribution. Fig. 24. To ascertain the extent to 


which the flow was uniform In Che aluminum test section, vertical surveys 
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taken along the side walls at axial locations of 77.5 and 83.31 inches* 
were compared ifith the vertical flow distribution at the same locations 
in the Mahogony test section. The uniformity of the Mach number distri- 
bution is far better with the aluminum test section. Therefore, the in- 
vestigations over the circular arc profile and the supercritical airfoil 
were conducted in the aluminum test section with the aluminum plates. 

The experimental results which were obtained are discussed below. 

6.1 Flow over Circular Arc Airfoil 

6.1.1 Model Solid Surface Test Results 

Before starting the investigation of the model porous surface 
effect on the flow field, pressure data were obtained with the no model 
porosity case and cosqpared to Llepmann's results, Ref. 16. Flguie 25 
shows the comparison between the two experimental results. The overall 
pressure coefficient distributions are similar with some differences due 
to the different airfoil chords used in the present study and in Llepman's 
work. Llepmann used a 3-lnch chord airfoil, whereas our model has a 
4-inch chord, both with 12-percent thickness. 

6.1.2 Preliminary Model Porous Surface 

A 12-percent-thlck circular arc airfoil with a 4-inch chord was 
placed in the center of the horizontal bottom wall of the 3-lnch by 15.4- 
inch wooden transonic test section for preliminary investigation. The 
model surface was made porous. Fig. 7b, with a cavity beneath, in the 
region extending from 60 to 90 percent of the chord length from the model 
leading edge. The results of the initial tests with the porous surface. 
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are presented in Figs. 26a and b, and indicate a change in the shock 
wave structure from a nearly normal shock wave to a Lambda shock wave, 
which is a combination of a normal and oblique shock waves. The oblique 
shock wave, caused by the flow circulation from the downstream to the up- 
stream of the shock wave through the cavity below, starts at the porous 
surface leading edge and extends upward to meet the normal shock wave at 
a point In the free stream above the porous surface. The change of the 
terminating shock wave into a normal one and the appearance of the new 
oblique shock wave at the porous surface leading edge show the effect of 
the porous surface and cavity on the flow field at transonic speeds. This 
effect, partially represented by the area between the oblique and normal 
shock waves. Increases with the free stream Mach number as shown in Figs. 
26a and b. These preliminary results were a positive sign and Indicated 
a need for further investigation of the concept of passive shock wave/ 
boundary layer interaction control. Because of its advantages discussed 
above, the aluminum test section was designed, constructed and used for 
the rest of the experiments. 

6.1.3 Effect of Aluminum Test Section Top Wall Porosity 

To study the effects of top wall porosity, Schlleren photographs 
of the flow field over the 12-percent- thick circular arc profile, mounted 
j.n the center of the first 20-lnch of the bottom wall, were taken for 
different wall porosities without the aluminum plates on the side walls 
of the test section. Figure 27a shows a strong, nearly normal shock wave 
terminating a large embedded supersonic region for the case of no top wall 
porosity. For the same control wedge setting at the test section exit. 
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Schlieren photograph of flow over 12-percent-thick 
circular arc airfoil, without top wall porosity 


Schlieren photograph of flow over 12-percent-thick 
circular arc airfoil, with 5 Z top wall porosity 
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an Incraaae In poroalcy decreaaas the size of the embedded supersonic 
zone and weakens the terminating shock wave and moves it upstream as 
indicated in Fig. 27b. As the porosity is increased to full top wall 
porosity of approximately 10-percent, the shock wave vanishes as shown 
in Fig. 27c. These preliminary tests in the aluminum test section were 
conducted with no model porosity. To minimize the effects of side wall 
deflection cn the flow field, the aluminum plates were added and the in- 
vestigation of the passive drag control of the circular arc airfoil was 
continued. 

6.1.4 Effect of Model Porosity With Large Cavity 
1.17 Percent Porosity 

The model surface was made porous by drilling .021-lnch diameter 
holes perpendicular to the surface. The full porous surface starts at 
60 percent of the chord from the model's leading edge and extends down- 
stream to 90 percent. The different porosities were developed through 
the chordwise variation of the porous surface length, by sealing or open- 
ing the appropriate holes, or by enlarging the hole size. This method 
enabled the investigation of the flow field with various porosities. The 
tests in this phase were conducted using a large cavity below the porous 
surface. The cavity was 2.125 inches long, 2.5 inches wide and 2 Inches 
deep. The 1.17 percent porous surface is located between 70 and 90 per- 
cent of the chord. Schlieren photographs and pressure distributions along 
the model centerline were obtained with and without the model porosity. 

The Schlieren photographs taken for the 1.17 percent porosity 
and the large cavity. Fig. 28b, show once again the effect of the model 
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Fig. 27c Schlieren photograph of flow over l.'?-percent- 
thlck circular arc airfoil, with 10 X top wall 
porosity 
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Fig. 28b Schlleren photograph of flow over 12-percent-thlck 
circular arc airfoil, with 1.17 Z porosity, M-0. 84 


Fig, 28a Schlleren photograph of flow over 12-percent-thlck 
circular arc airfoil , without porosity, M-0.84 
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porosity on chs shock wavs scructurs. An obllqus shock wavs, appaarlng 
at chs porous surface's leading edge, owacs the normal shock wavs in a 
Lambda configuration. The initial shock wavs without porosity. Fig. 28a, 
has bean affected and its shape changed into a normal shock wave. Fig. 

28b. However, the location of uie terminating shock wave remains un> 
changed. The shadowgraph photograph. Fig. 28c, taken with the 1.17 per- 
cent model porosity, shows the same shock wave configuration. 

The model Mach number distribution without the porosity and j 

with 1.17 percent porosity indicates that the porosity decreases the | 

local Mach number over the portion of the porous surface upstream of the 

shock wave location, while increasing the Mach number downstream. j 

1.6 Percent Porosity | 

The initial porosity of 1.17 percent was brought to 1.6 percent 
by increasing the length of the model porous surface, starting at 57 per' 

i 

cent of Che chord and ending at 80 percent. Experiments were conducted 
with and without porosity and the large cavity. Schlleren photographs, I 

pressure distributions over the model, and wake surveys, were obtained. i 

The Schlleren photograph of the flow field over the 12-percent-thlck cir- 
cular arc airfoil with and without porosity shorn the same results as in 
the case of 1.17 percent porosity. Because of the porosity effect the 
shock wave moved slightly forward with a slight decrease in height. A 
related phenomena obtained from the model pressure measurements shows a 
decrease of the maximim local Mach number and a shift in the shock wave 
location, with porosity. The Mach number distribution remains imchanged 
upstream and downstream of the porous surface. A slight recovery of Mach 

j 

J 

] 
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number Is noticed over the portion of the porous surface downstream of 
the shock wave. For a lower free stream Mach number, the surface Mach 
number distribution Is less affected by the porosity than at higher 
speeds. However, there Is still a slight decrease of the maximum local 
Mach number with porosity. 

The wake survey data presented In Fig. 29, were obtained by 
measuring the Impact pressures within the wake with a rake. The rake 
centerline Is located at a distance of approximately 40 percent of the 
chord downstream of the model's trailing edge, and has 10 Impact pressure 
probes. The wake survey starts at the test section bottom surface and 
extends upward to around 1.75 Inches. The total pressure ratio distribu- 
tion, shown In Fig. 29, Is obtained by measuring the Impact pressure at 
different vertical locations In the wake and dividing It by the total 
pressure of the potential flow field measured in the settling chamber. 
Figure 29 shows the effect of the model porosity on the impact pressure 
distribution. The losses of the total pressures across the shock wave are 
less with porosity than without, except within the boundary layer between 
0 and .30 inch from the bottom surface. 

To investigate the effects of model porosity on the free stream 
Mach number, vertical surveys of the Mach number distributions at two 
different locations ahead of the model were obtained. The vertical Mach 
number distributions were determined from the static pressures measured 
on the test section side wall at different vertical locations. The side 
wall pressure data with and without porosity, shows that there Is no effect 
of the model porosity on the free stream Mach number. 
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0.94 Percent Porosltv 

The porosity was made smaller by decreasing the length of the 
model's porous surface so that It started it 61 percent of the chord from 
the leading edge and ended at 71 percent with the same hole slse. The 
Investigation of the passive drag control was conducted with the large 
cavity to study the effects of different porosities on the flow field 
without changing the cavity size beneath the porous surface. Schlleren 
photographs obtained without and with 0.94 percent model porosity, show 
the decrease In the shock wave height due to the effect of the porosity. 
The same shape Lambda shock wave was expected, but the brightness of the 
photograph made It impossible to confirm. 

0, .94 and 1.6 Percent Porosity 

Three of the previous porosity cases are compared for the same 
free scream Mach number and the result of the model Mach number distribu- 
tions Is presented in Fig. 30. As discussed previously, the no porosity 
case gives a higher maximum local Mach number. However, the case of 1.6 
percent porosity has less effect on Che Mach number distribution chan the 
0.94 percent case for the same cavity size. Upstream and downstream of 
Che porous surface, the distribution Is unchanged for all three cases. 

For Che three porosity cases, the differences In the Mach number distribu- 
tions aroimd the mid-chord get smaller as the free scream Mach number de- 
creases . 


6.1.5 Effect of Cavity Length 

Keeping Che same 2-lnch depth, Che cavity was made smaller by 


decreasing Its length In the chordwlse direction. The length was changed 









from 2.125 Inches Co 1.69 Inches, making Che new cavicy dimensions 1.69 
inches long by 2.50 inches wide by 2 inches deep. The scudy of che cavicy 
size is mainly for praccical reasons. Since che space wichin che wLsag 
of an airplane i^ere che caviCy will be insCalled is limiced, a smaller 
cavicy is preferable. Experimencal daca from che previous large cavicy 
and che new smaller one are compared for che same model poroslcy size 
of 0.94 percenC and che same free scream Mach number. The resulcs are 
presenced in Figs. 31a and b. Figure 31a i,hows Che model Mach number 
dlscrlbuclons for a conCrol wedge deflecclon ac Che Cesc secclon exic of 
0 degrees and 0.94 percenc model poroslcy. Boch dlscrlbuclons are similar 
ups cream and downscream of Che porous surface, buc over Ic and around 
Che model mid-chord, che smaller cavicy gives higher dlsCribuCion which 
corresponds Co a higher life. Figure 31b shows che impacc pressure dls- 
crlbucion measured in Che wake downscream of Che model Cralllng edge. 
Compared co Che large cavlCy case, che losses of impacc pressure Chrough 
Che shock wave is sllghcly less wlCh che smaller cavicy. Based on Chase 
resulcs, che small cavicy was selecced for furcher InvesclgaClon of Che 
passive drag concrol concepC. 

6.1.6 EffecC of Cavicy Depch 

The lengch and widch of che small cavicy are kepc Che same, in 
Che following series of experimenCs, buc Che depch was changed from 2 
Inches Co 0.25 inch. Schlieren pho cographs, pressure dlscrlbuclons over 
Che model, and wake surveys were obcalned wich 0.94 percenC porosiCy Co 
scudy Che ef feces of che cavicy depch on Che flow field. Figure 32 shows 
Che Schlieren pho Cograph of che flow field over che 12-percenC- Chick 
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circular arc. This time the Lambda shock wave shape is quite visible. 

The shock wave terminating the supersonic region is a straight normal 
shock wave of approximately 2-lnch height. The oblique shock wave rising 
from the porous surface leading edge extends upward and meets the ter> 
mlnatlng shock wave at a distance of 0.70 inch from the model surface. 

The two black marks shown in the picture locate the chord of the airfoil, 
one mark at the model's trailing edge and one, not quite visible, at the 
leading edge. 

The effect of the cavity depth on the model Mach number distri> 
butlon is presented in Fig. 33a with a 0.94 percent model porosity and 
a 5 degree wedge deflection at the test section exit. The 2-inch depth 
produces a lower Mach number distribution which corresponds to a lift 
loss. This is a disadvantage in the present investigation because the 
objective is to achieve drag reduction with minimum lift losses. Figure 
33b shows no difference in the Impact pressure distribution from the wake 
survey for both depth cases. The surveys were made with 0.94 percent 
model porosity and zero degree wedge deflection at the test section exit. 
Based on these results, the 0.25 inch depth was selected since it mini- 
mizes the lift losses with the porous surface model. 

6.1.7 Effect of Large Model Porosity with Small Cavity 

The 0.94 percent model porosity was increased to 2.04 percent 
by enlarging the diameter of the porous surface holes from 0.021 inch to 
0.031 inch. The pressure distribution over the model as well as the im- 
pact pressures from the wake were obtained. The results were compared 
against the no porosity and the previous 0.94 percent porosity case. 


n 
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Fig. 33a Mach number distributions over 12-percent 
thick circular arc airfoil with 0.25 and 2 
inch cavity depth, 0.94 percent porosity, 
M-0.83 


Fig. 33b Wake total pressure distributions for 12-percent 
thick circular arc airfoil with 0.25 and 2-lnch 
cavity depth, 0.94 percent porosity, t(,*0.83 






Figure 34a shows Che pressure coefficient dlscrlbucion over the model 
for 0.824 free stream Mach number. From the model's leading edge up to 
30 percent of the chord* the pressure distribution remains unaffected by 
Che porosity. Between 30 to 60 percent of the chord, both 0.94 percent 
and 2.04 percent porosities slightly Increase the pressure coefficient 
valxies similarly. Beyond that point, the 3 distributions align again 
with slightly lower values for the 2.04 percent case. The tests were con> 
ducted with the removal of the boundary layer so that the new boundary 
layer would start at the airfoil leading edge stagnation point. A 0.22 
inch thick rake with 10 Impact pressure probes, was used for the wake 
survey. Figure 34b shows the point drag distribution for the three model 
porosities of 0, 0.94 and 2.04 percent. The point drag coefficient is 
derived using Eq. (26). For a free stream Mach number of 0.827, the 
porous model surface Induces a slight drag reduction; however, increasing 
the porosity from 0.94 to 2.04 percent has negligible effects on the drag 
distribution, as shown in Fig. 34b. 

6.1.8 Effect of the Rake Blockage 

Because of the blockage effects at transonic speeds, where most 
of the actual tests are conducted, the thick rake was replaced by a thin 
one so that the effect of the rake would be minimized. The tests indi- 
cated the effect of the rake size on the model Mach number distribution. 
The thinner rake yielded the higher model Mach number diatrlbutlon, al- 
though the distribution remained unaffected over the first 30 percent of 
the chord. The tests were conducted with a free stream Mach number of 
0.828, a model porosity of 2.04 percent, and the small cavity of 0.25 




Fig. 34a Pressure coefficient distributions over 12~percent* 
thick circular arc airfoil without .with 0.94 and 2.04 
percent porosity .and l/4>inch deep cevlty. t^O.83 


mur OM amxttm vnoo 

Fig. 34b Local drag coefficient distributions for 12-perccnt 
thick circular arc airfoil without .with 0.94 and 2.04 
percent porosity, and 1/4-inch deep cavity 
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Inch deep. 


6,1.9 Effect of thn Thin Rake Position 

The thin rake selected for the reaalnlng experinental tests 
consists of 8 Impact pressure probes of .031 Inch diameter each, trlth a 
total height of .875 Inch and a thickness of .063 Inch. For the complete 
experimental data within the wake downstream of the model, the vertical 
survey was made In 2 steps. The rake was first positioned at the test 
section bottom surface and then was moved vertically to a position such 
that It covered the necessary height for a complete w^e survey, over* 
lapping a couple of the impact pressure probes between the 2 positions. 

To check the effect of the rake's vertical position on the flow field, 
pressure data were :«btalned with a model porosity of 2.04 percent and 
with the removal of the boundary layer. The tests showed practically no 
effect of the rake position on the model Mach number distribution. 

6.2 Flow Over Supercritical Airfoil 

12 

A 14-percent- thick NASA supercritical airfoil, with a 4-lnch 
chord was constructed and used In the Investigation of passive drag con- 
trol. The supercritical airfoils, now being used for the current gen- 
eration of aircraft, are shaped to reduce the drag associated with energy 
losses due to the shock wave and flow separation at transonic speeds. 
However, despite this advantage of drag reduction over the conventional 
airfoils, a strong viscous Interaction between the shock wave and the 
boundary layer still occurs, and the need to control It Is siore necessary 
today than It has been In the past to achieve a new generation of more 
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efficient alrcrefc. The research described In this section will focus 
on this problem. The present supercritical airfoil, designed for a 
normal force coefficient of 0.7, was mounted on the middle of the first 
20-lnch test section bottom surface. Seventeen pressure taps, spresd 
along the model surface centerline, were installed to messure the static 
pressure distribution. Before the model surface was made porous, pre- 
liminary tests were conducted with the smooth surface. Schlleren photo- 
graphs and pressure distributions, as well as wake surveys, were obtained 
for a range of free stream Mach numbers. Similar tests, with the excep- 
tion of the Schlleren photographs, were conducted by NASA on the sane 
type of airfoil with the smooth surface corresponding to the no porosity 
case. The comparison of the results is discussed in the following sec- 
tion. 
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6.2.1 Smooth Model Surface Test Results 

In the tests conducted in the 3- Inch by 15.4-lnch RPI test sec- 
tion over the 14-percent-thick supercritical alr'lell at a free stream Mach 
number of 0.84, the boundary layer removal has no effect on the shock 
wave location, as evidenced by the Schlleren photographs. Figure 35 shows 
expansion waves in the supersonic region over the model surface, terminated 
by a strong shock wave of approximately 2.46 inch in height, and located 
at 77 percent of the chore. The expansion waves, departing from the air- 
foil surface, meet the sonic line bounding the embedded supersonic sons 
and are reflected as compression waves back to the model surface. The 
expansion waves are di.verging, while the compression waves are converging. 
When these compression waves finally coalesce, they form a strong shock 
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Fig. 35 Schlieren photograph of flow over 14-percent- 
thlck smooth surface supercritical alrfoi' , 
>fe0.804 
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wave which terminates the supersonic region. This strong shock wave of 
large height and the expansion waves starting within 23 percent of the 
chord from the leading edge, define a large embedded supersonic region, 
bounded by the sonic line, the model surface, and the terminating shock 
wave. 

Static pressures were measured on the surface of the model and 

used to determine the local Mach number distribution for a range of free 

stream Mach numbers. The surface pressure measurements were obtained 

from a chordwlse row of orifices located along the model centerline. 

Orifices were concentrated near the airfoil leading edge to define the 

stagnation point as well as the pressure gradient In this region. 

The chordwlse Mach number distribution, shown In Fig. 36, Indl- 

18 

cates the similarity between NASA's results and the present experimental 
data, with a descrepancy In the similarity around the trailing edge, where 
the NASA values of the Mach numbers are lower. Both Mach number distri- 
butions show the same maxlmtim local Mach number of 1.26 at the same loca- 
tion, and the same shock wave location. NASA experiments In this case 
were conducted with a free stream Mach number of 0.78 and a normal force 

coefficient C of 0.51. However, with a lower free stream Mach number of 
n 

0.74 and a C of 0.26, the two Mach number distributions are shown to be 
n 

more alike, even In the trailing edge region. Fig. 37. Here, the maximum 
local Mach number Is unity, and the graph of the distributions Is flatter. 
Indicating a more constant distribution than In the previous case of M^ 
equal to 0.78. The model was mounted on the test section bottom wall. 

Fig. 8b. Consequently, only the upper surface and the leading edge were 


Fig. 36 Mach number distributions over R.P.I and NASA 
I4-percent-thlck supercritical airfoils without 
porosity, t|*0.80 
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Fig. 37 Mach number distributions over R.P.I and NASA 
14-percent-thlck supercritical airfoils without 
porosity, M-0.74 
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investigated, with only one pressure measurement taken on the model 
bottom surface at a distance of 3 percent of the chord from the leading 
edge. The Mach number at this point was compared to the NASA bottom sur- 
face distribution for a free scream Mach nuisd>er of 0.74. Because of the 

deficiency in the removal of the boundary layer at the model leading edge, 

18 

the presently obtained value is higher chan NASA's. The removal was 
not adequate; therefore, Che new boundary layer did not start at the 
model's stagnation point. Fig. 38 . 

As Che free stream Mach number decreases to .65, Fig. 39 indi- 
cates a more constant Mach number value of .81, from near the leading 
edge to 70 percent of the chord. Compared to NASA's results, both dis- 
tributions coincide practically from the leading edge to Che trailing 
edge. Again, for a free stream Mach number of 0.65, the value of the 
pressure coefficient obtained on the airfoil bottom surface near Che 
leading edge is lower compared to Che NASA bottom surface results for 
the same reasons of the boundary layer removal as discussed previously. 

Impact pressure data within Che wake downstream of the model's 
trailing edge were obtained for different free stream Mach numbers using 
the Chin rake. To obtain the section drag coefficient, Eq. (25), the 
point drag coefficients were computed for each total pressure measure- 
ment in the wake, by using Eq. (26). These point drag coefficients were 
Chen summed by numerical integration across the wake, based on the trape- 
zoidal method. Figure 40a shows a good drag distribution over a wide 
range of free scream Mach numbers, with no measurable shock losses up to 
Che drag divergence Mach number of approximately 0.76. Figure 40a also 
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shows Che NASA results wich a slightly higher drag divergence Mach number. 
The variation of the section drag coefficient with the free stream Mach 
number, with and without model porosity is presented In Fig. 40b. 

The overall Mach number distributions for a range of free scream 
Mach numbers show In general good agreement wich the results presented in 
Fig. 18. From these experimental results, the location of Che porous sur- 
face on the 14-percent- thick supercritical airfoil was selected so that 
Che shock wave would be standing In the middle. The experiments on the 
previous 4-lnch chord circular arc airfoil were used as a guide for this 
next phase of the supercritical airfoil Investigation. Consequently, 2.5 
percent porosity with a 3/4-lnch cavity beneath It \taa selected, as well 
as Che thin rake for the wake survey downstream of the airfoil trailing 
edge. The experimental results are discussed In Che following section. 

6.2.2 Effect of Model Porosity with Large Cavi ty 
2.5 Percent Porosity 

The model's surface was made porous by drilling 0.025 inch dia- 
meter holes perpendicular to the surface. The porous surface, consisting 
of 18 spanwlse rows with 38 holes each, is located between 56 percent and 
S3 percent of the chord from the airfoil leading edge. Based on the model 
surface area, the porosity is defined as the ratio between the total area 
of Che porous surface holes and the area of the model. Consequently, Che 
corresponding porosity of the above geometry is 2.5 percent. The cavity 
beneath the porous surface was made 3/4 inch deep, 3 Inches wide and 1.08 
Inches long. The lengths of the cavity and the porous surface, in the 
chordwise direction, are the same for the case of 2.5 percent porosity. 

The holes of the porous surface can be either sealed or enlarged to obtain 




the desired porosity* The model's surface is made non-porous by sealing 
all the holes. Schlieren photographs, pressure distributions and wake 
surveys were obtained without and with the 2.5 percent porosity for a 
range of free stream Mach numbers. 

The Schlieren photographs taken with a free stream Hach number 
of .806, shown in Figs. 41a and b, indicate a sizable change in the shock 
wave structure due to the effect of the model's porosity. As discussed 
previously, the porous surface with the cavity permits a part of the bound- 
ary layer to move from downstream to the upstream of the shock wave, send- 
ing compression waves which produce an oblique shock wave. The oblique 
shock wave originating from the porous surface leading edge reaches the 
terminating shock wave at a distance of 1.23 Inches from the model sur- 
face. The location of the terminating shock wave in the middle of the 
porous surface, 80 percent from the leading edge, was not affected by 
the porosity; its shape, however, was changed into a normal shock wave. 
Expansion waves can be seen in both Schlieren photographs, with and with- 
out the porosity, upstream of the shock wave near the model leading edge, 
indicating the presence of a large embedded supersonic region. The thick 
black mark in the pictures indicates the location of the model's trailing 
edge. The leading edge is not quite visible because of the stress con- 
centration in that region of the plexiglass side walls. The two thin 
black marks approximate the location of the porous surface on the model. 

As the free stream Mach number decreases to 0.78, the Schlieren 
photographs. Figs. 42a and b, without and with the 2.5 percent porosity 
respectively, show the effect of the porosity on the location of the shock 
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Schlieren photograph of flow over 14-percent-thick 
supercritical airfoil, without porosity, 806 


Schlieren photograph of flow over 14-percent-thick 
supercritical airfoil, with 2.5Z porosity, 3/4-lnch 
cavlty,M«0.806 
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Fig. 42a Schlieren photograph of flow over 14-percent-thick 
supercritical airfoil, without porosity ,M-0. 78 


Fig. 42b Schlieren photograph of flow over 14-percent-thlci 
supercritical airfoil, with 2.5Z porosity, 3/4-lnch 
cavity, M;0.78 






wave and its strength. The shock wave moved upstream from 62 to 57 per- 
cent of the chord. Its shape, however, remained a single normal shock 
wave with a decrease In the height from 1.23 to 1.10 Inches, corresponding 
to a reduction In Its strength, as well as a reduction In the else of the 
embedded supersonic region. The shock wave was still over the porous sur- 
face. With further reductions In the free stream Mach number, the shock 
wave disappeared, as shown In Fig. 43 for a Mach number of 0.74. 

The Mach number distributions were obtained from the pressure 
measurements over the model surface. The distributions with 2.5 percent 
porosity and without porosity, at, respectively, 0.806 and .803 free stream 
Mach numbers, were ccmpared. The Mach number distribution remained un- 
affected by the porosity up to the mid-chord location. The porosity 
effects began to appear when the values of the Mach number were decreased 
up to the location of the shock wave at nearly 80 percent of the chord. 

The maximum Mach nuisber, located close to the mid -chord, was decreased 
from 1.23 to 1.18. Downstream of the shock wave, the Mach number was 
higher with the porosity than without, which corresponded to a decrease 
In the local pressure, and led to a reduced flow separation. For better 
comparison the Mach number distributions of 2.5 percent and no porosity 
cases were compared again, but this time for the same free stream Mach 
number of 0.806, as shown In Fig. 44a. The same observations as before 
were made here, except that the no porosity case gave a higher maximiun 
Mach number of 1.26 Instead of 1.23 due to a higher M^, and that the two 
distributions aligned at the niodel's trailing edge. The corresponding 
pressure coefficient distributions are shown In Fig. 44b. 
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Fig. 43 Schlieren photograph of over 14-percent-thick. 

supercritical airfoil without porosity ,M«0. 74 



Flg.A4a Mach number diatrlbutlona over lA>percenc-thlck 
supercritical airfoil without and with 2.5 percent 
porosity, 3/4-inch deep cavity. 1^0.806 


Fig. 44b Pressure coefficient distributions over 14-percent 
thick supercritical airfoil without and with 2.5 
percent porosity, 3/4-inch deep cavity, 1^0.806 
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For A Mach number of 0.78, the reaulta praaancad In Fig. 45a 
ahow alnllar affacta of the poroua aurfaca on the Mach number diacrlbu- 
cion, aa wall aa for the higher free stream Mach number of 0.806. The 
dlatrlbuclon remained unaffected from the leading edge to the model &.id- 
chord, and Chen decreaaed over the porous surface upatream of the shock 
wave location at around 60 percent of the chord. Downstream of the shock 
wave, Che porosity Increased the Mach number and islnlmlzed the flow aeparao 
tlon phenomenon. 

For a Mach number of 0.74, the effects of the porosity on the 
Mach number distribution started 25 percent farther upatream than the case 
of higher Mach number and extended downstream to the end of the porous 
surface. The effects were then insignificant until uhe trailing edge, 
trtiere the Mach number became higher with porosity, as shown In Fig. 45b. 

The flow was fully subscnlc over the model surface, with a relatively 
constant Mach nuoiber distribution between 20 and 50 percent of the chord. 

As the free scream Mach number was decreased to 0*65, the Mach 
nusd>er distributions, shown In Fig. 45c, became flatter between 17 and 
60 percent of the chord, with a nearly constant value of 0.81 for the no 
porosity case, and 0.79 for the 2.5 percent case. The coi^arlson between 
the two distributions remained the same as In Che previous higher Mach 
number condition. 

To investigate the effect of the airfoil porous surface on the 
drag reduction, drag coefficient distributions were obtained from impact 
pressure steasurements In the wake downstream of the model trailing edge 
as shown in Figs. 4 and 9. The experiments were conducted for a series 
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Fig* 45a Mach number distributions over 14-percent- 
thick supercritical airfoil without and with 
2*5 percent porosity, 3/4-inch cavity ,^0*78 
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Fig* 45b Mach number distributions over 14-percent- 
thick supercritical airfoil without and with 
2.5 percent porosity, 3 /4-inch cavity ,^0.74 
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Fig. 45c Mach number distributions over 14-percent- 
thick supercritical airfoil without and with 
2.5 percent porosity, 3/4-inch deep cavity. 
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of free stream Mach numbers. For a free stream Mach number of 0.807, 

25 

the point drag coefficient distributions. Indicated a lower drag dis- 
tribution with 2.5 percent porosity than with the no porosity case through- 
out the survey height, except within the boundary layer between 0 and 
0.25 inch. The difference in the local drag of the porous and non-porous 
surfaces decreased as the height Increases. The two drag distributions 
finally align at a height of approximately 1.75 inches. The numerical 
integration of the point drag coefficients across the 1.75 inch height 
showed a net drag reduction due to the airfoil porous surface. The tests, 
repeated for approximately the same free stream Mach number of 0.806, 
showed similar drag coefficient distributions, with a drag reduction. The 
same observation was obtained, from Fig. 46a, with the free stream Mach 
number of 0.803 for the no porosity case and 0.806 for the 2.5 percent 
porosity one. Even though the Mach number with the no porosity case was 
lower than with the 2.5 percent case, there was still a net drag reduc- 
tion over the airfoil upper surface of 15 percent due to the model poros- 
ity. The corresponding impact pressure distributions are shown in Fig. 

46b. 

As the free stream Mach number was decreased to 0.78 for both 
cases, with and without porosity, there was a drag reduction with porosity, 
except within the first 0.25 inch of the wake survey height where the drag 
was higher. Fig. 46c . The point drag distribution decreased as the height 
Increases and finally vanished at a height of 0.80 inch for both with and 
without porosity cases. At this free stream Mach number of 0.78, the net 
drag was higher with the porous surface. 
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Ac a Mach number of 0.741, the results presented in Fig. 46d 
i 'Hlcated a higher drag with the porous surface. The two point drag co- 
efficient distributions decreased rapidly within Che first 0.25 Inch and 
vanish for approximately the rest of the height, indicating a shock-free 
flow over the model as shown in Che corresponding Schlleren photograph. 
Fig. 43, discussed earlier. 

Consequently, at transonic Mach nimbers Inducing a strong shock 
wave, the 2.5 percent model porous surface with a 3/4-lnch cavity does 
reduce the profile drag over the airfoil. However, at subsonic speeds 
at which the shock wave vanishes and where there Is no wave drag, the drag 
Is higher. 

More data with the 3/4-lnch deep cavity were obtained, but for 
a smaller porous surface of 1.25 percent, the results are presented in 
Che following section. 

1.25 Percent Porosity 

To Investigate the effect of porosity size on Che shock wave/ 
boundary layer interaction and on the flow field, the 2.5 percent surface 
porosity was reduced to 1.25 percent by sealing every ocher hole. The 
cavity dimensions were kept the same. Experimental tests were conducted 
for a series of free stream Mach numbers ranging from transonic to sub- 
sonic speeds. Schlleren photographs, pressure distributions and wake 
surveys were obtained with the new porosity and same cavity size. 

The Schlleren photograph taken for M^ of 0.807, shows a normal 
shock wave with a 2.4 luch height, terminating the local supersonic region 
at 80 percent of the chord. A series of compression waves emanating from 
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Che porous surface mec Che cemlnaclng shock wave wichln a heighc of 1.22 
inches. Expansion waves were also seen wichin Che supersonic zone down- 
stream of 24 percenc of Che chord. Their expansion into Che flow field 
as well as Che heighC of Che cerminaclng shock wave Indlcaced Che presence 
of a large local supersonic region. Compared Co Che no porosicy case. 

Fig. 47a, Che change of Che shock wave shape inco a normal one and Che 
appearance of Che compression waves over che porous surface, indlcaCed 
Che effecC of Che porosicy on Che flow field. However, che effecc on Che 
cerminaclng shock wave locaclon was negligible. 

For a sllghcly lower free scream Mach number of 0.804, che 
effecc of che compression waves origins cing from che porous surface on 
che shape of Che Cerminaclng shock wave was nociceable, as shown in Figs. 
47a and b. Beside Che change inco a normal shape, Che compression waves 
had a tendency Co pull upstream che point of Incersecclon with Che cer- 
minaclng shock wave. Expansion waves close Co Che model leading edge and 
a normal Cerminaclng shock wave of 2.34 Inches heighc located at 78 per- 
cenc of che chord indicated che presence of a large local supersonic 
region embedded in a subsonic flow. The effecc of che porous surface on 
che cerminaclng shock wave location was negligible, as indicated in Figs. 
47a and b. 

For a subsonic free stream Mach number of 0.78, Che porosicy 
effecc on che terminating shock wave locaclon, at Che beglimlng of che 
porous surface, was negligible. However, a decrease in che heighc from 
1.43 to 1.17 inches was produced as shown in Figs. 47c and 42a. Expansion 
waves between che 37 percent chord point and Che normal shock wave defined 
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Fig. 47a Schlleren photograph of flow over 14-percent-thick 

supercritical air foil, without porosity, 3/4-lnch cavity, 
M- 0.803 

oo 


Fig. 47b Schlleren photograph of flow over 14-percent-thlck 
supercritical airfoil, with 1.252 porosity, 3/4-lnch 
cavity ,M“0. 804 








a small local supersonic region as indicated in Fig. 47c. Hence, for a 
lower free stream Mach number of 0.78, the shock wave became weaker under 
the porosity effect, but did not take on the Lambda configuration. Also, 
its location remained approximately the same. At 0 Mach number of 0.743, 
the flow was entirely subsonic over the model without any shock wave. 

The data of the model Mach number distribution for different 
free stream Mach numbers were obtained. Fig. 48a shows the distributions 
for M^ equal to .803 and .807, without and with 1.25 percent porosity, 
respectively. The distribution remained unaffected by the porosity from 
the model leading edge to 65 percent of the chord. In this region the 
Mach number increased over the model surface up to a maximum value of 
1.23 near the mid-chord. The sudden decrease in the Mach number through 
the shock wave, located at approximately 78 percent of the chord, was 
less using the porous surface, which corresponds to a lower static pres- 
sure Jump through the shock wave. Referring to Eq. (15), the ratio of 
the statin pressures before and after the shock wave defined its strength. 
Since the Mach number ahead of the model is the same with and without 
porosity, a lower static pressure ratio indicates a weaker shock wave. 

The same conclusion can be drawn from Eq. (16). The static pressure 
difference p2*’Pi through the shock wave was less with the porosity than 
without. For the same static pressure, p^, the corresponding increase 
in the entropy will be less according to Eq. (17). This means that the 
porous surface induced a weaker shock wave over the airfoil. 

The pressure distributions data for the cases of no porosity 
and of 1.25 percent porosity with a 3/4-inch cavity, were obtained for 



128 



I 

I 


I 

I 

i 

1 

I 

[ 

r. 

[ 



i: 


free scream Mach numbers of 0.803 and 0.804, respectively. The resulcs 
In terms of pressure coefficients are presented in Fig. 48b. The dis- 
tribution over Che model remained unaffected up to the mid-chord point; 
then Che porous surface slightly Increased the pressure coefficient between 
approximately the model mid-chord and the shock wave location at 78 per- 
cent of the chord. Downstream of the shock wave to the trailing edge, 
the porosity Induced a slightly lower pressure coefficient, C^, compensa- 
ting for the gain in upstream of the shock wave. 

18 

The drag data for the 14-percent-thlck supercritical airfoil, 
with 1.25 percent porosity and without porosity, at respectively M^ of 
0.807 and 0.803, are presented in Fig. 48c. Notice that the point drag 
coefficient distributions decreased rapidly within the first 0.50 inch 
height, with slightly higher values for the 1.25 percent porosity. How- 
ever, between 0.50 inch and the rest of the survey height, the local drag 
coefficient was lower, corresponding to a reduction in the wave drag. Al- 
though the free stream Mach number was higher with porosity than without, 
the integration of the point drag coefficients showed a slight net drag | 

reduction with the porous surface. 

To investigate the porosity size effect on the shock wave, pres- 
sure distribution and the wake intact pressures, the two previous cases 
of different porosities were compared for the same 3/4-lnch cavity size. 

The results are discussed below. 

2.5 Percent and 1.25» Percent Porosity 

1 

The data for the full porosity case corresponding to 2.5 per- j 

1 

i 

cent and of the 1.25 percent case were compared for the same cavity size | 
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Fig. 48a Mach number distributions over 14-percent~thick 

supercritical airfoil without and with 1.2S percent 
porosity, 3/4-inch deep cavity 
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Fig. 48b Pressure coefficient distributions over 14-percent 
thick supercritical airfoil without and with 1.25 
percent porosity, 3/4-lnch deep cavity 
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Local drag coafficiene dlatrlbutloos for 
14-percant-chlck aupercrltical airfoil 
tflehout and with 1.25 percent poroaity* 
3/4-lncb deep cavity 
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of 3/4»lach depth. The two Hech nuaber dletrlbutlons for of 0.806 
presented in Fig. A9e» Indicate that the 2.5 percent porosity has loiMr 
Mach nuaber values over the porous surface and downstream* except for the 
area becveen 70 and 80 percent of the chord, where the 1.25 percent case 
was slightly lower. Consequently, the increase in porosity decreased 
the local maximum Mach number. 

The point drag coefficient distributions from the wake survey, 
shown in Fig. 49b for a free atream Mach number of 0.806, indicated a 
lower drag with the 2.5 percent over the 1.25 percent porosity throughout 
the wake survey height. 

The experimental reaults of the no porosity case, as well as 
the 2.5 percent and 1.25 percent porosity cases with 3/4<-lnch cavity, are 
compared for different free stream Mach numbers in the following section. 
Zero Percent. 2.5 Percent and 1.25 Percent Porosity 

The data concerning the Mach number distribution and drag were 
compared between the three porosity cases of zero percent, 2.5 percent 
and 1.25 percent porosity for a range of free stream Mach numbers. 

Figure 50a shows the Mach number distribution obtained with M^ 
of 0.803, 0.804 and 0.806 for zero, 2.5 and 1.25 percent porosity, re- 
spectively. The porosity did not affect the distribution over the first 
half of the model surface. Between 50 percent and approximately 70 per- 
cent, corresponding to the porous region upstream of the shock wave, the 
porosity decreased Che Mach number; the higher the porosity, the lower 
the Mach number. However, downstream of the shock wave location, between 
78 percent of the chord sxji the model trailing edge, the porosity slightly 
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Increased the Mach number; the higher the porosity, the higher the Mach 
number. The Mach number distributions showed vlrtiially no effect of the 
porosity on the shock wave location. 

The Mach number distributions for the three previous porosity 
cases were obtained with the same free stream Mach number of 0.78. As 
shown In Fig. 50b, the distribution remained unaffected up to approximately 
the beginning of the model porous surface which Is located between 56 and 
83 percent of the chord. Over the section of the porous surface upstream 
of the shock wave, the Increase In porosity decreased the Mach number. 
However, over the section downstream from the shock wave, the increase in 
porosity Increased the Mach number. Downstream from the porous surface, 
the distributions of the three porosity cases coincided, although the 2.5 
percent case did have a slightly higher value. The shock wave moved up- 
stream from approximately 62 percent of the chord with no porosity to 58 
percent with porosity. The same observation was obtained from the corre- 
sponding Schlleren photographs. Consequently, the partial movement of 
the decelerated flow from the back of the shock wave to the front Increased 
the Mach number downstream from the shock wave, and sent compression waves 
ahead of the shock wave, which decreased the upstream Mach number. As 
the porosity Increased the amount of the displaced flow was larger, thereby 
Increasing Its effect on altering the values of the Mach numbers, as seen 
In Fig. 50b. 

At M • .74, the flow over the model was entirely subsonic, as 
shown in Fig. 50c. The Mach nianber distribution remained unaffected by 
the porosity for the first 50 percent of the chord. Between approximately 
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50 and 70 percent of the chord corresponding to the upstream portion of 
the porous surface, the Increase In porosity decreased the Mach number. 
However, over the rest of the porous area, the Mach number was slightly 
higher with 2.5 percent porosity, and remained unaffected with 1.25 per- 
cent porosity. Downstream from the porous surface, over the region lo- 
cated between 87 percent of the chord and the airfoil trailing edge, the 
Increase In porosity Increased the Mach number. Even In the absence of 
the shock wave, the flow over the downstream portion of the porous surface 
carried a higher pressure than the flow over the upstream part. Therefore, 
a circulation of part of the flow within the boundary layer was established 
between the downstream and the upstream portions of the porous surface. 

This circulation through the porous surface, which moved a decelerated 
flow from one side and Injected It Into the boundary layer ahead of the 
shock wave. Induced respectively an Increase and a decrease in the existing 
Mach numbers over these regions. By this process, the boundary layer sepa- 
ration was minimized, leading to an Increase In the Mach number near the 
trailing edge, as shown In Fig. 50c. 

The experimental data for the drag, obtained from the Impact 
pressure wake surveys for the three porosity cases of 0, 1.25 and 2.5 per- 
cent, are presented in Fig. 51a. The tests were conducted with the 3/4- 
Inch cavity and free stream Mach numbers of 0.803, 0.804 and 0.806 respec- 
tively. Within Che first 0.25 Inch of Che boundary layer, the porosity 
led to a slightly higher local drag coefficient, with no difference between 
the 1.25 and 2.5 percent cases. However, for Che rest of the wake survey 
height, the Increase In porosity decreased the local drag coefficient. 
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The reduction of the wave drag was greater with the 2.5 percent porosity, 
eventhough the corresponding free stream Mach number was slightly higher 
than the 1.25 percent and the no porosity cases. The three drag dlstrl- 
butlons coincided at a height of approximately 1.65 Inches. The numerical 
Integration of the point drag coefficients throughout the wake survey 
height, showed a net drag reduction over the airfoil upper surface of 
approximately 14 percent with the 2.5 percent porosity, and 5 percent 
reduction for the 1.25 percent case. Corresponding total pressure ratio 
distributions are presented In Fig. 51b. 

For a lower free stream Mach number of 0.78, Fig. 51c, the In- 
crease In porosity increased the local drag coefficient within the first 
0.25 Inch of the wake survey height. However, over the remaining height, 
the 1.25 percent porosity led to a better drag reduction than the 2.5 per- 
cent case. For this low free stream Mach number, the porosity Induced a 
reduction In the wave drag, but a slight increase In the total drag. 

A further decrease In the free stream Mach number to 0.74 es- 
tablished an entire subsonic flow over the model without wave drag. Within 
the boundary layer, the increase In porosity Increased the local drag, as 
shown In figs. 51d. A corresponding total pressure ratio increase Is 
shown In Fig. 51e. 

6.2.3 Effect of Model Porosity with Small Cavity 

For further investigation of the passive drag control concept, 
the depth of the cavity beneath the porous surface was reduced and ex- 
perimental data with the 2.5 percent and 1.25 percent porosity were ob- 
tained. The results will be discussed in the next two sections. 
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Flg.Sld Local drag coefficient distributions for 14-percent* 
thick supercritical airfoil without, with 1.25 and 
2.5 percent porosity, 3/4-lnch deep cavity 


Fig.Sle Wake total pressure distributions for 14-percent- 
thick supercritical airfoil without, with 1.25 and 
2.5 percent porosity, 3/4-lnch deep cavity 
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2.5 Ptrcent Porosltv 

To Investigate the effect of the cavity size on the shock wave/ 
boundary layer Interaction and on the flow field, the 3/4-lnch cavity 
depth was reduced to 1/4-lnch. The rest of the cavity dimensions re» 
malned the same. The 2.5 percent porosity was selected first, because 
It Induces a higher drag reduction at high Mach numbers. 

Schlleren photographs and pressure data, as well as Impact 
pressure wake surveys, were obtained for a range of free stream Mach 
numbers. The Schlleren photograph, (-.g. 52a, obtained with M^ equal to 
0.807, shows the effect of the porosity and cavity, on the shock wave 
shape. 

The difference In pressure across the terminating shock wave 
caused a part of the decelerated flow downstream from It to flow upstream 
through the porous surface, sending compression waves Into that region. 
Consequently, a new oblique shock wave was produced at the porous surface 
leading edge, \dilch joins the terminating shock wave at a point 0.73 Inch 
above the model surface. The effect of the porosity and the cavity on 
the shock wave terminating the supersonic region was visible compared to 
the no porosity case presented in Fig. 41b. The shock wave with porosity 
was normal and weak over the 0.73-lnch portion close to the model surface. 
Under the porosity effect, the shock wave moved from 79 percent of the 
chord to 73 percent. Expansion waves could be seen between the 24 percent 
chord point and the normal terminating shock wave. Directly downstream 
from this normal shock wave, at a location of 78 percent of the chord 
from the leading edge, a very weak shock wave of approximately .58 Inch 
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Fig. 52a Schlieren photograph of flow over 14-percent-thick 
supercritical airfoil, with 2,571 porosity, 1/4-inch 
deep cavity,M=0.807 


Fig. 52b Schlieren photograph of flow over 14-percent-thlck 
supercritical airfoil. with 2.5% porosity, 1/4-lnch 
deep cavity, 1^=0.803 
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occurred over the porous surface. The Lanbda shape seeaed to be a char* 
accerlsclc of the porosity effect on the shock wave at transonic speeds. 

The tests were repeated under the same conditions with equal 
to 0.803, Fig. 52b. This figure shows a series of strong con^resslon 
waves upstream of the terminating shock wave, produced by the circulation 
of the removed flow through the porous surface. The effect of those com- 
pression waves on the flow field became obvious by slaq>ly observing the 
perfectly normal termlnai-^ng shock wave. Conq>ared to the no porosity 
case. Fig. 47a, the shock wave moved slightly upstream to a location of 
approximately 77.6 percent of the chord Instead of 79 percent, while the 
shock wave height of 2.42 Inches remained the same. Expansion waves within 
the embedded supersonic region are visible In the Schlleren photograph. 

With equal to 0.78, Fig. S2c shows a weak shock wave at 
approximately S3 percent of the chord, whereas with no model porosity, 
the location was at 62 percent of the chord. Fig. 42a. The shock wave 
height was decreased from 1.45 Inches to 0.73 Inch, as was the size of 
the expansion waves. Consequently, the size of the supersonic region, 
as well as the local supersonic Mach numbers, were reduced under the 
porosity effect. 

The Schlleren photograph obtained for a subsonic free stream 
Mach number of 0.74 Indicates an entirely subsonic flow over the model 
without any shock wave. The boundary layer Is not visible because the 
photograph was taken with a vertical knife edge. Attempts have been made 
(though without much success) to obtain Schlleren photographs with a 
r horizontal knife edge to show the boundary layer behavior. The deficiency 
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Fig. 52c Schlieren photograph of flow over 14-percent- 
thlck supercritical airfoil, 2. 5Z porosity, 
1/4-lnch cavity, M-C 78 
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waa cauaad nalnly by cha thlckneaa of cha aide wall plaxlglaaa which 
creacea a diacortlon of che incldaot light bean through the teat aectlon. 
Similar reaulta trare obtained with a free atream Mach number < f 0.6S. 

Preaaure data were obtained for a range of free atre a Mach 
numbera. Fig. S3a ahowa the conparlaon of the model Mach number dlatrl- 
butlon without and with 2.5 percent poroalty for reapectlve free atream 
Mach aumbera of 0.806 and 0.807. From the model'a leading edge to the 
mid-chord, the diatrlbutlon remained unaffected. However, from SO percent 
of the chord to the trailing edge, the Mach number waa lower In the caae 
uf 2.5 percent porosity. The effect of the porosity was greater over 
the porous surface. The difference In the Mach number acrosa the shock 
wave was smaller with the 2.5 percent porosity, corresponding to a lower 
Jump In pressure across the shock wave. Since the Mach number, or the 
corresponding pressure, Just downstream of the shock wave was the same 
for both with and without porosity, as Indicated in Fig. 53a, the ratio 
of the static pressure, shock wave was lower ifith 2.5 

percent porosity. The static pressure ratio, 

the strength of the shock wave, and thus the shock wave was weaker with 
the porosity. A weaker shock wave implied a lower entropy Increase, as 
Indicated by Eq. (17), or a lower total pressure loss, Eq. (18), causing 
a wave drag reduction. The corresponding distributions of the pressure 
coefficient data are presented In Fig. 53b. A similar observation Is 
obtained with a free stream Mach number of 0.804. 

For a lower M of 0.78, Fig. 53c, the effect of the porosity 
on the decrease of the Mach number was less severe than at a higher free 






stream Mach number. The Mach number distribution with porosity, was 
slightly lower than without porosity between approximately 33 percent 
of the chord and the shock wave location at 53 percent and slightly higher 
between S3 percent and 65 percent of the chord. From 65 percent of the 
chord to the trailing edge, the porous case was again lower, with approxi- 
mately a constant difference between the two cases. 

At a free stream Mach number of 0.74, the flow was entirely sub- 
sonic, as shown In Fig. 53d, and the effect of the porosity extended to 
the leading edge. However, the effect was greater over the porous surface 
location. The two distributions with and without porosity coincided be- 
tween approximately 77 percent of the chord and 90 percent. At the trail- 
ing edge the porosity Induced a higher Mach number. A similar observation 
of the porosity effect on the Mach number distribution was obtained at a 
free stream Mach number of 0.65. However, the porosity In this case was 
slightly less effective. For a transonic Mach number of 0.807, the data 
on the drag obtained from the impact pressure wake survey shows a drag 
reduction with 2.5 percent porosity throughout the wake survey height. 

The local drag coefficient vanished at approximately 1.30 Inches height 
with porosity, while with no porosity. It was still non-zero at that 
location. Consequently, the strength of the shock wave as well as Its 
height were decreased by the porosity. The Integration of the local drag 
coefficients with and without porosity Indicated a net upper surface drag 
reduction of 35 percent. 

For a free stream Mach number of 0.804, the 2.5 percent porosity 
and the 1/4-lnch cavity reduced the drag throughout the wake irvey 






height, as shown in Fig. S4a. The drag reduction was minimum within the 
first 0.25 Inch above the model surface, and a maximum between 0.25 Inch 
and 0.50 inch, and then decreased within the rest of the wake survey height. 
The local drag coefficient with porosity vanished at 1.70 Inches from the 
airfoil surface, a shorter distance than without porosity, corresponding 
to a reduction In the height of the shock wave. The summation of the 
local drag coefficients throughout the wake height showed a total upper 
surface drag reduction of approximately 38 percent. A corresponding re- 
duction In total pressure due to the porous surface Is shown In Fig. 54b. 

The decrease In the free stream Mach number to 0.78 resulted In 
higher drag with 2.5 percent porosity than without porosity over the 0.25 
Inch region Just above the model surface, as shown In Fig. 54c. However, 
for the rest of the height, the porosity reduced the drag. The point 
drag coefficient vanished at 1.00 Inch height In the no porosity case and 
at 0.40 inch with the 2.5 percent case. The calculation Indicates a slight 
Increase in the total drag, although there Is a local wave drag reduction. 

For a lower cf 0.74, the drag was higher with the porosity 
up to a height of 0.30 Inch from the airfoil surface, as shown In Fig. 54d. 
Within this height, the local drag coefficient for both with and without 
porosity, decreased and vanished beyond that height limit. A similar drag 
distribution was obtained with a free stream Mach number of 0.65. 

The results obtained so far demonstrated that the concept of the 
passive drag reduction at transonic speeds Is truly promising, but seeking 
a better porous surface effectiveness, further Investigation and testing 
were done with a 1.25 percent porosity. 
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Flg.SAa Local drag coefficient distributions for 14-percent 
thick supercritical airfoil without and with 2.5 
percent porosity, 1/4-lnch deep cavity 
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Fig. 54b Wake total pressure distrlbutloas for 14-percent 
thick supercritical airfoil without and with 2.5 
percent porosity,!/ 4-inch deep cavity 
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Fig. 54c Local drag coefficient distributions for 14>percent 
thick supercritical airfoil without and with 2.5 
percent porosity, 1/4-inch deep cavity 
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Flg.54d Local drag coefficient distributions for 14-percent 
thick supercritical airfoil without and with 2.5 
percent porosity, 1/4-inch deep cavity 



1.25 Percent Porosity 


The 2.5 percent porosity was reduced to 1.25 percent by sealing 
every other hole of the model porous surface. The tests were conducted 
with the 1/4-lnch cavity depth. Schlleren photographs were obtained for 
a range of free stream Mach numbers. The effect of the porosity on the 
shock wave can be seen by referring to the Schlleren photograph In Fig. 

55a, obtained for a Mach number 0.804. Compression waves produced by the 
flow circulation through the porous surface emerge upstream of the ter- 
minating shock wave and affected Its shape. Expansion curves between the 
airfoil leading edge and the shock wave were visible In the saaie figure. 

The shock wave height Indicated the presence of a large local supersonic 
region over the airfoil surface. The 1.11-lnch portion of the shock wave 
close to the model surface became normal under the compression waves effect. 
Consequently, the entire terminating shock wave changed from a nearly 
normal shape, as seen In Fig. 4 /a for the no porosity case, to a completely 
normal one. The location of the shock wave at 79 percent of the chord, as 
well as its height, remained unaffected by the porosity. 

For a lower subsonic Mach number of 0.78, the Schlleren photo- 
graph in Fig. 55b showed a small supersonic region defined by the expansion 
waves and a small normal shock wave. The height of this terminating shock 
wave decreased from 1.45 Inches to 0.97 inch by the porosity effect, and 
Its location moved slightly forward from 63 percent of the chord to 61 
percent. 

The Schlleren photographs obtained tor free stream Mach numbers 
of 0.74 and 0.65 show the flow field over the modal surface without any 
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Fig. 55a Schlieren photograph of flow over 14-percent-thick 
supercritical airfoil, with 1.25 % porosity , 1/4-inch 
deep acvity, I’3*0.304 
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Fig, 53b Schlieren photograph cf flow over 14-percent--thick 

supercritical airfoil, with 1.25 % porosity, 1/4-inch 
deep cavity,M=»0,78 
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shock wave. The flow la entirely subsonic and the boundary layer behavior 
Is the only optical Information left which can predict the porosity effect. 
Unfortunately, this boundary layer Is not visible on the Schlleren photo- 
graphs taken with a vertical knife edge. Attempts have been made to obtain 
It with the knife edge In the horizontal position. The thickness of the 
side wall plexiglass produces a non-uniform contrast of the flow field 
picture with the horizontal knife edge. Therefore, the Schlleren photo- 
graphs obtained did not provide enough Information about the boundary layer. 
Improvements on this situation will be made In the near future. 

6.2.4 Effect of Cavity Depth 

To Investigate the most effective of the two cavity depths, Hach 
number distributions as well as local drag coefficient distributions were 
compaZvid for the 3/4-lnch and 1/4-lnch cavities with the same model porosity 
of 2.5 and 1.25 percent. 

2.5 Percent Porosity 

The Mach number distributions for M of 0.806 and 0.807 are pre- 

OP 

sented In Fig. 56a. The distributions were similar upstream of the porous 
surface. However, from the porous surface leading edge to the model trail- 
ing edge, the 1/4-lnch depth Induced a lower Mach number. The difference 
In the Mach number was larger at the trailing edge. The corresponding 
drag distributions for the same free stream conditions are presented In 
Fig, 56b. The 1/4-lnch cavity depth gave a lower drag than the 3/4-lnch. 

As the vertical height Increased, the drag decreased to a zero value at a 
height of 1.30 Inches for the 1/4-lnch case, while It was still non-zero 
for the 3/4-lnch case. 


Fig. 56a Mach number distributions over 14-percent-thlck 
supercritical airfoil with 3/4-lnch and 1/4-inch 
deep cavity, 2. 5 percent porosity 
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Fig. 56b Local drag coefficient distributions for 14-percent 
thick supercritical airfoil vrlth 3/4-inch and 1/4- 
inch deep cavity, 2. 5 percent porosity 











157 


For a free screaa Mach nunbar of 0.78, Che data of Che drag 
presenced In Fig. 56c, indicace Chat Che l/4»lnch cavlcy Induced a lover 
drag dlscrlbuclon chan Che 3/4-lnch depch. The local drag coefficient of 
Che l/4>lnch depch vanished at 0.40 Inch while chac of Che 3/4-lnch van- 
ished ac 0.7S Inch. 

The results showed chat there was a better drag reduction with 
Che 1/4-lnch cavity depth, with little loss of lift, as compared to the 
3/4-lnch case. To verify this deduction, both depth cases were compared 
again, but this time for a different model porosity of 1.25 percent. 

1.25 Percent Porosity 

The comparison of the Mach number distributions obtained with 
of 0.807 Is presented in Fig. 57a. The distributions td.ch the 3/4-lnch 
and the 1/4-lnch cavity were similar over the first hslf of the airfoil 
profile as well as over the last 20 percent of the chord. Between 50 and 
65 percent of the chord, the 1/4-lnch cavity caused a lower Mach number 
than the 3/4-lnch. However, between 65 and 80 percent. It Induced a 
higher Mach number. Indeed, the loss In the corresponding lift obtained 
between 50 and 65 percent of the chord was compensated by the gain In 
lift between 65 and 80 percent. 

For a free stream Mach number of 0.78, Fig. 57b shows absolutely 
no difference between the Mach number distribution for these two cavity 
depths. Comparison of the drag distribution between the two cavity depth 
for M^ of 0.804 is presented in Fig. S7c. The 1/4-lnch cavity produced 
a lower drag throughout the wake survey height. A similar result was 
obtained with the free stream Mach numbers of 0.78, Fig. 57d, and 0.74, 




Fig. 57a Mach number distributions for I4«percent- 
thick supercritical airfoil with 3/4-lnch 
and 1/4-lnch deep cavity, 1.25 percent porosity 




Fig. 57b Mach number distributions over 14-percent- 
thlck supercritical airfoil with 3/4-lnch 
and 1/4-lnch deep cavity, 1.25 percent porosity 





Fig. 57c Local dra; coaf^iclmt dlatributlona for 14- 
parcent-chlck supercritical airfoil with 3/4- 
inch and 1/4-lnch deep cavity, 1.25 7 , poroalty 


Flg.57d Local drag coefficient dlatributlona for 14 
percent supercritical airfoil with 3/4-lnch 
and 1/4-lnch deep cavity ,1.25% porosity 
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Flg.57e Local drag coefficient distributions for 
14-percent-thlck supercritical airfoil 
with 3/4-inch and I /4-inch deep cavity, 
1.25 percent porosity 
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Fig. S7e. With the 1.25 percent porosity, the 1/4-lnch deep cavity pro- 
duced a lower drag with no loss In lift as compared to the 3/4-lnch case. 

6.2.5 Effect of Humidity 

To Investigate the humidity effect on the flow field, the dryer 
located at the entrance of the Transonic Wind Tunnel, as shown In Fig. 1, 
was reactivated to bring the relative humidity from around 50 percent down 
to around 23 percent. Schlleren photographs, pressure distributions and 
Impact pressure wake survey obtained without model porosity, before and 
after the dryer reactivation were compared. 

To best simulate free flight conditions In the Transonic Wind 
Tunnel, the humidity was brought to a lower range by reactivating the 
dryer for a period of approximately 48 hours. The Schlleren photograph 
In Fig. 58a, obtained for a free stream Mach number of 0.806, shows a 
displaced shock wave from 79 percent to 76 percent of the chord due to 
the decrease in the relative humidity. However, the shock wave height 
and shape remained unaffected. 

The Mach number distributions before and after dryer reactiva- 
tion are compared In Fig. 58b, with no model porosity and with a free 
stream Mach number of 0.806. The decrease In humidity increased the 
local maximum Mach number located at approximately 55 percent of the chord, 
and produced a slightly lower Mach number distribution between 65 percent 
of the chord and the airfoil trailing edge. However, over the first 
ualf of the model chord, the distribution was unaffected. 

The corresponding drag distributions obtained under the same 
conditions are presented In Fig. 58c. The Increase In humidity decreased 
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Fig. 58b Mach number distributions over 14-percent- 

thick supercritical airfoil before and after 
dryer reactivation, without porosity ,M>0. 806 
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Fig. 58c Local drag coefficient distributions' for 14- 

percent-thlck supercritical airfoil before and 
after dryer reactivation, without porosity ,F^0. 806 
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the drag throughout the height of the wake survey and brought It to zero 
at a smaller height of 1.45 Inches. The difference In the drag between 
the two distributions before and after the dryer reactivation was neg‘> 
llglble within the first 0.25 Inch, and thereafter Increased with In- 
creasing height. 

To Investigate the humidity effect on the drag reduction, data 
were obtained after the dryer reactivation, with 2.5 percent porosity and 
1/4-lnch deep cavity as well as with no porosity. Tho choice of the 2.5 
percent porosity and 1/4-lnch cavity sizes for a comparison with the no 
porosity case was based on the previous test results. 

A Schlleren photograph obtained with the 2.5 percent porosity 
and a free stream Mach number 0.806 is shown In Fig. 59a. Compared to 
the no porosity case. Fig. 58a, the terminating shock wave was once again 
affected by the porosity In its shape as well as In Its height. New cont- 
presslon waves produced by the flow circulation through the porous surface 
as discussed previously, transformed the existing shock wave into a normal 
one. A decrease in the shock wave height from 2.28 Inches to 2.13 Inches 
was also produced. However, the shock wave location at 76 percent of 
the chord remained unaffected by the porosity. The compression waves 
departing from the model porous surface reached the terminating shock 
wave within a height of 0.97 Inch. 

The Mach number distribution under the same conditions as above 
Is presented In Fig. 59b. The distribution remained unaffected by the 


2.5 percent porosity from the model leading edge to approximately half 
of *the chord. Between 50 and 70 percent of the chord, the porosity de- 






Fig « 59b Mach number distributions over 14-percent* 

thick supercritical airfoil without and with 
2.5 percent porosity,! /4-inch cavity after 
dryer reactivation 


Fig. 59c Pressure coefficient distributions over 14- 
percent-thlck supercritical airfoil without 
and with 2.5 Z porosity, 1/4-lnch cavity, after 
dryer reactivation 





167 


creased the Mach cumber, vhlle between 70 and 90 percent, it increased 
it. The two distributions of with and without porosity coincided at the 
airfoil trailing edge. The porosity decreased the Mach number over the 
portion of the porous surface upstream of the shock wave and Increased 
it downstream. The corresponding pressure coefficient distributions are 
presented in Fig. 5 9c. 

The drag distribution with the previous test conditions, is 
presented in Fig. 59d. The 2.5 percent porosity caused a lower drag dis- 
tribution throughout the survey height, with approximately zero drag be- 
yond 0.50 inch height. The maximum drag reduction occurred between the 
heights of approximately 0.25 and 0.75 inch. This reduction was mainly 
a wave drag reduction, which occurred within the compression waves in- 
fluence region. The summation of the local drag coefficients sho%ied a 
net upper surface drag reduction of approximately 18 percent. The 
corresponding total pressure distribution is shown in Fig. 59e. 

6.2.6 Effect of Test Section Top Wall Porosity 

The preceding tests were conducted with the test section porous 
top wall fully closed. To investigate the effect of this top wall porosity 
on the flow field and the drag reduction, half of the full porosity was 
selected. The full porosity, based on the total area of the top wall, 
was approximately 10 percent. 

Schlleren photograph, pressure distributions and Impact pressure 
wake surveys were obtained without and with 2.5% porosity and 1/4-lnch ca- 
vity for a free stream Mach number 0.806. The Schlleren photograph 
obtained with no model porosity is presented in Fig. 60a. A strong shock 
wave located at 76 percent of the chord, with a 2.23-inch height termlxiated 
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Flg.59d Local drag coefficient distributions for 
lA-percent-tblck supercritical airfoil 
without and with 2.5 percent porosity, 1/4 
inch deep cavity, after dryer reactivation 


Fig..S9e Wake total pressure distributions for 14- 

percent-thlck supercritical airfoil without 
and with 2.5 X porosity, 1 /4-inch deep cavity 
after dryer reactivation 
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60s Schlteren photograph of flow over 14-percent-thlck 
supercritical airfoil .without porosity, 1/2 top wall 
porosity ,MgO. 806 


Fig. 60b Schlieren photograph of flow over 14-percent-thlck 

supercritical airfoil, 2. 5 Xporosity and 1/4-lnch cavity 
1/2 top wall porosity, t^O. 806 
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a large supersonic zone within which expansion waves could be seen. The 
effect of the porosity Is evident In Fig. 60b, where the shock wa\e has 
been changed Into a normal one under the Influence of the compression 
waves. The shock wave location was not affected by either the model or 
the top wall poror.^ui . However, its height decreased to 2 inches due 
to the model poro&lty effect. 

With the same half top wall porosity and of .806, the Mach 
number distributions were obtained with and without model porosity. The 
trends shown In Fig. 60c are similar to those In the previous Hacb number 
figures. The distribution remained unaffected by the model porosity over 
the first half of the airfoil surface. Between 50 and 70 percent of the 
chord, the model porosity decreased the Mach number while between 70 per- 
cent and the airfoil trailing edge. It slightly Increased it. The de- 
crease In the Mach number was produced over the portion of the porous 
surface upstream of the shock wave, while the increase was downstream. 

The corresponding pressure coefficient distribution Is presented In Fig. 
60d. 

The drag distributions with 2.5 percent and without model poro- 
sity for a free stream Mach number of .806 are presented In Fig. 60e. 

The porosity Induced a drag reduction thrcu'/hout the wake survey height. 
The reduction was again maximum between approximately 0.25 and 0.75 Inch, 
the region of the compression waves' influence. The estimated net upper 
surface drag reduction was 8 percent. 

6.2.7 Analysis of Shock Wave Losses 

To confirm the experimental wave drag reduction results, a 


IJ 


II 


Fig. 60c Mach number distributions over 14-percent- 

thick supercritical airfoil without and with 
2.5 X porosity. 1/4-inch cavity, 1/2 top wall 
porosity, M«0.806 
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Fig.60d Pressure coefficient distributions over 14- 
percent-thick supercritical airfoil without 
and with 2.5 % porosity, 1/4-inch cavity, 

1/2 top wall porosity , IJ^O.806 
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Fig.60e Local drag coefflclant dlstrlbuciona for 14- 
pcrceat-chick supercritical airfoil without 
and with 2.5 percent porosity, 1/4-lnch deep 
cavity, 1/2 top wall porosity, {j[*0.806 
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brief theoretical aDelysls wee applied to the flow acroaa the shock wave 

over the model surface with and without porosity. By assuming. In the 

case of no porosity, that the terminating shock wave has a normal shi^e, 

the normal shock wave relatione indicate the following results 
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where and P^^ total pressures ahead and behind the shock wave; 
P^ and ^2 static pressures and and S 2 the entropy per unit mass. 

With the porous surface model at transonic Mach number, the ex> 
periments showed the following Lambda Shock Wave System consisting of an 
oblique and a normal shock wave. 


Oblique shock 
wave 


...Normal shock wave 
Pprous surface 
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Using the oblique shock wave relations between 1 and 1', and the normal 
relations between 1* and 2, we obtain 
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where M! Is given by Eq. (20. 
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If we assume the flow deflection angle 6 after the shock wave to be neg- 
ligible compared to the shock wave angle B« we obtain 
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The entropy increase is given by the following equations 
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Consequently, the entropy increase between 1 and 2 is given by 
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The static pressure increase is given by 
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This yields the pressure Increase across the shock wave system 


(15) 
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By using the experimental value of the local Mach number 
ahead of the shock wave, for the cases of with and without porosity, and 
by measuring B from the Schlleren photographs, the above equations give 
the following results: 


Porosity Porositi 


^ 02^^01 

P2/P1 


.9857 


1.686 


.9989 

.0011 

1.402 


It follows from these results that the total pressure loss, the entropy 
Increase and the static pressure Increase are all less with porosity than 
without. Indeed, the theory confirms the experimental wave drag reduc- 
tion results. The flow circulation through the porous surface and the 
cavity had provided a new oblique shock wave which lowered the local Mach 
number Just upstream of the shock wave terminating the embedded super- 
sonic region, thus weakening the strength of the shock. Figs. 61a and b 
with and without porosity respectively show a good agreement between the 
theoretical and experimental total pressure loss results. 

The 1/4-lnch deep cavity was expected to provide larger pressure 


difference between the flow upstream of the terminating shock wave and 
the flow within the cavity than the 3/4-lnch cavity. This was confirmed 
by the cavity pressure measurements. Therefore, the smaller cavity produces 
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a stronger oblique shock wave and consequently, weakens the normal shock 
wave strength more than the larger cavity. Since the main total pressure 
losses are produced across the normal shock wave rather than across the 
oblique shock, the weaker terminating shock wave associated with the 1/4- 
Inch cavity produces larger drag reduction, as shown In Fig. 56b. 

Although there Is an undesirable compression effect on the lift 
due to the oblique shock wave upstream of the terminating shock wave, the 
porous surface has, on the other hand, the desirable lift Increase effect 
over the region downstream from the shock, caused by the decelerated flow 
circulation as shown In Fig. S9b. This flow mechanism i^ch produces a 
lift decrease upstream of the shock wave and an Increase downstream from 
It tends to keep the net lift unchanged with and without porosity. 

For smaller model porosity of 1.25 percent, compression waves 
are produced Instead of an oblique shock wave, lowering the local Mach 
number Just ahead of the terminating shock, but not as low as with the 
oblique shock In the 2.5 percent porosity case. Therefore, the drag 
reduction with the 1.25 percent model porosity Is less, as shown In Fig. 
51a. 

The data with Che 2.5 percent porosity and 1/4-lnch cavity showed 
an appreciable drag reduction associated with Che Lambda Shock Wave System 
and with the decelerated flow circulation with or without a negligible 
loss in lift. The section drag distribution with and without porosity, 
presented In Fig. 40b, shows the porosity effect at transonic Mach numbers, 
causing either a drag reduction or a speed increase. 



PART 7 


CONaUSION 

A 3-lnch by 15.4-lnch Transonic Wind Tunnel was designed, 
constructed, and calibrated to operate over a Mach number range of 0.5 
to 0.95. Modifications of the test section were made to account for the 
boimdary layer growth along the tunnel walls. The flow Mach number In 
the test section was controlled by the adjustable wedges located on the 
top and bottom walls of the diffuser. 

The bottom wall boundary layer just ahead of the circular 
arc or supercritical airfoil was removed by a transverse slot across 
the wall. A valve was Installed in the suction piping to regulate the 
amount of boundary layer that was removed. 

To minimize the wall Interference effects on the flow field 
over the airfoil, a variable porosity top wall was installed. By sli- 
ding the movable backside of two porous plates while the alrslde plate 
remained stationary. It was possible to vary the opening from zero to 
10 percent. 

A quick-acting pneumatic valve, located downstream from the 
diffuser section, was used to establish the flow in approximately one 
second. The static and Impact pressures In the test section and over 
the models were measured by mercury manometers, which had a solenoid 
valve for each manometer bank. The duration of the test flow, and the 
recording of the static pressures, and the Schlleren photographs were 
all controlled by a timer system. Thus, the pressures and the Schlleren 
photograph were taken simultaneously. 
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Luclte side walls were Installed on the test section, over 
idilch were fastened thick aluminum plates with cut-outs to permit the 
observation of the boundary layer and the shock waves over the airfoils 
by the Schlleren system. The Schlleren photographs were taken by a ca- 
mera which was triggered by the timer system used to operate the tunnel. 

It was possible to take the Schlleren photographs with either a horlson- 
tal or a vertical knife edge position. 

The impact pressure dlstrlbutlone downstream from the airfoil 
trailing edge were measured with a narrow cross-sectional rake of 8 im- 
pact preasue probes. With this rake It was possible to measure the ver- 
tical Impact pressures and calculate the upper surface profile drag for 
the various porosities, cavity depths, and flow Mach numbers. 

Investigations of the passive drag control concept were con- 
ducted with a 12-percent-thlck circular arc and a 14-percent-thlck su- 
percritical airfoil, which were placed on the bottom wall of the test 
section. The porous surface with the cavity beneath it was positioned 
at the shock wave location at transonic Mach numbers. The porosities 
of 1.25 and 2.5 percent, over a large and a small cavity, were investi- 
gated to determine the optimum size of porosity and cavity. 

At transonic speeds, the airfoil porous surface with the 
cavity beneath It produced a lambda shock wave system consisting of an 
oblique and a normal shock wave. The oblique shock wave lowered the lo- 
cal Mach number just upstream the normal shock terminating the emLedded 
supersonic region, thus weakening the strength of the normal shock wave. 

In addition, the circulation of the decelerated flow through the porous 
surface from the downstream to the upstream terminating shock wave location 
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lowered the pressure distribution over the downstream region, minimi- 
zing the flow separation. 

The wake Impact pressure surveys showed an appreciably lo- 
wer impact pressure loss and a drag reduction for the airfoils with 
the porous surface, with or without a negligible loss in the lift. 

A higher drag reduction was obtained with 2.S percent porosity and the 
small cavity at transonic Mach numbers. At low subsonic speeds where 
there was no shock wave over the airfoil surface, the porosity increa- 
sed the drag. To overcome this shortcoming at low Mach numbers, an 
anticipated solution Is to seal the porous surface. 

To determine the overall merit of the concept and to apply 
the idea in future aircraft generations, more development and testing 
must be done. A full-scale program la also needed to demonstrate that 
the concept Is truly ready for civil and military applications. Such 
full-scale programs usually involve flight testing, which Is beyond our 
apabllitles. 

Many economic factors may limit technological development, 
but the present concept appears to be economically feasible because of 
Its simplicity, and the fact that It needs no additional power to func- 
tion. 

The results of the present research show that such a passive 
drag reduction concept will be feasible for applications' to future air- 
craft design, since any drag reduction represents either a direct saving 
in fuel or an Increase in performance. 
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